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Abstract 

This  project  endeavors  to  find  whether  it  is  feasible  to  use  an  increase  in  surface 
area  as  a  way  of  increasing  the  drag  on  an  orbiting  object,  thus  decreasing  its  orbital 
lifetime.  The  surface  area  increase  can  be  achieved  by  an  apparatus  that  deploys  a 
balloon.  The  balloon  will  act  as  a  parachute  that  will  decrease  the  potential  energy 
of  the  object  through  atmospheric  drag.  This  is  most  effective  by  objects  that  reach 
the  Low  Earth  altitudes  of  less  than  500  kilometers,  where  an  object  is  encountering 
a  firmer  atmospheric  density. 

The  project  is  carried  out  through  propagating  three  different  element  sets  to 
reentry  using  STK®.  The  orbital  paths  generated  by  the  software  are  then  graphed 
in  Microsoft  Excel®  and  presented. 

The  analysis  is  divided  into  four  main  studies.  The  first  study  focuses  on  con¬ 
firming  the  effects  the  atmospheric  instability  has  on  the  long  term  predictions  of  a 
natural  decay.  The  second  study  explores  how  an  increase  in  the  scale  of  the  drag, 
at  different  points  in  the  orbital  path,  affects  the  reentry  time.  The  third  study  in¬ 
vestigates  a  specific  increase  of  area  to  mass  ratio  (A/M)  at  different  points  in  the 
trajectory.  This  is  to  survey  changes  in  the  variability  of  the  reentry  time  prediction 
and  how  the  reentry  location’s  variability  is  altered.  To  finish  off,  the  last  study 
examines  how  A/M  manipulation  affects  the  reentry  time  prediction. 

The  project  discovered  that  an  increase  in  A/M  decreases  the  variability  in 
reentry  prediction.  Furthermore,  it  discerns  the  exponential  relationship  between  the 
time  to  reentry  and  the  A/M. 
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I.  Introduction 

Since  the  launch  of  Sputnik  I  the  space  community  has  placed  every  possible  object 
into  orbit,  from  animals  to  sophisticated  satellites.  These  activities  have  had 
important  missions  but  there  has  not  been  a  strong  and  consistent  worldwide  effort 
into  mitigating  the  space  debris  problem  we  are  facing  today.  Understandably,  the 
space  debris  problem  was  left  alone  with  the  purpose  of  benefiting  the  advancement 
of  technology,  as  well  as  not  adding  another  burden  to  the  space  acquisitions  and 
budgets  of  the  world.  Nonetheless,  space  operational  regions  are  a  limited  resource 
that  must  be  protected. 

1.1  Space  Environment 

Currently  the  total  debris  amount  hovers  about  the  twelve  thousand  mark, 
which  includes  objects  in  all  orbits.  Projections  (using  current  trends  and  not  ac¬ 
counting  for  mitigation  measures)  for  Low  Earth  Orbiting  debris  only  are  in  the 
eleven  thousand  range  by  the  year  2010.  To  make  things  worse,  some  satellites’  mis¬ 
sions  never  even  get  started.  In  the  last  few  decades,  total  satellite  failures  have 
happened  with  an  almost  predictable  regularity.  These  failures  have  left  the  objects 
floating  in  orbit  to  naturally  decay  back  to  the  earth.  Some  of  these  failures  have  been 
of  objects  close  to  earth  with  periods  small  enough  to  reenter  the  earth’s  atmosphere 
within  a  few  years  or  even  months. 

1.2  Removal  of  Space  Debris  Needed 

Objects  in  orbit  that  have  no  control  system  are,  in  essence,  kinetic  weapons 
with  hypervelocities  in  the  range  of  kilometers  per  second.  For  this  reason,  it  is 


1 


advantageous  to  have  the  objects  reenter  since  it  lessens  the  possibility  of  orbital 
collisions  amongst  debris  and  decreases  the  danger  of  damage  or  destruction  of  active 
satellites.  The  root  concept  behind  this  thesis  is  the  basic  principal  that  drag  causes 
an  object  to  decrease  in  orbital  period,  eventually  making  it  reenter. 

Another  reason  to  have  an  object  reenter  in  a  controlled  manner  is  the  safety 
of  the  population.  There  are  satellites  that  carry  large  payloads  such  as  the  Skylab 
reentry  on  July  11,  1979,  which  was  around  one  hundred  tons.  There  are  dangerous 
satellites  such  as  Cosmos  954  that  reentered  on  January  24,  1978,  which  had  a  nuclear 
reactor  and  spread  radioactive  fuel  on  a  600  kilometer  path.  Other  failed  payloads 
include  USA  193,  which  was  shattered  to  pieces  by  a  missile  on  February  21,  2008,  in 
the  interest  of  population  safety. 

1.3  Thesis  Purpose 

The  goal  of  this  thesis  is  to  bring  forth  and  establish  the  idea  that  we  can  gain 
a  certain  amount  of  control  over  reentering  objects  that  have  no  control.  The  idea 
would  be  implemented  by  a  small  independent  system  that  will  only  manipulate  the 
drag  of  the  object.  Consequently,  the  time  and  place  of  reentry  within  a  particular 
window  may  be  chosen.  Fundamentally,  this  thesis  presents  the  proposal  of  increasing 
the  atmospheric  drag  of  an  object  through  the  increase  in  the  object’s  surface  area, 
thus  increasing  its  area-to-mass  ratio.  The  surface  area  increase  will  be  through  the 
inflation  of  a  balloon  that  will  act  as  a  parachute  on  the  object.  The  new  area-to-mass 
ratio  is  what  will  directly  increase  the  atmospheric  drag  on  the  object,  decreasing  its 
period,  and  ultimately  speeding  its  reentry. 

1.4  Atmospheric  Instability 

One  huge  problem  of  reentry  prediction  is  the  amount  of  instability  in  the  earth’s 
atmosphere.  This  instability  makes  the  atmosphere  difficult  to  model,  making  any 
type  of  prediction  that  is  further  than  two  days  into  the  future  extremely  unreliable. 
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In  effect,  another  thing  gained  from  decreasing  the  time  the  object  is  left  in  orbit  is 
that  we  can  make  the  reentry  prediction  much  more  reliable. 

1.5  Previous  Attempts  at  Enforced  Reentry 

Intentional  reentry  caused  by  an  increase  in  drag  has  been  proposed  before,  but 
through  using  tethers  on  satellites.  Such  proposal  did  not  include  gaining  control 
over  the  time  and  place  of  reentry.  The  use  of  a  tether  has  very  much  the  same  effect 
as  the  proposed  parachute  with  the  distinct  difference  that  a  parachute  will  increase 
the  radius  of  the  object’s  cross-sectional  area  by,  at  most,  tens  of  meters.  A  tether, 
on  the  other  hand,  will  typically  have  the  length  of  several  kilometers  and  a  width 
of  millimeters.  These  last  details  increase  the  possibility  of  something  going  wrong 
during  the  orbital  lifetime  of  the  object. 
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II.  Problem  Background  and  Review 


This  section  will  address  the  increasing  problem  of  debris  build,  the  rising  hazard 
to  the  population  from  debris  reentering  the  earth,  and  how  the  unpredictability 
of  the  atmosphere  makes  it  difficult  for  us  to  deal  with  reentry  prediction.  It  will  also 
mention  how  balloon  satellites  were  already  used  in  the  1960s.  To  conclude,  it  will 
bring  up  the  “square-cube”  law  may  be  used  to  help  with  the  space  debris  problem. 

2.1  A  Space  Debris  Problem 

For  the  last  few  decades  the  United  States  (US),  as  well  as  other  nations,  has 
launched  into  space  satellites,  experiments,  and  even  parts  for  a  laboratory,  i.e.  the 
International  Space  Station,  leaving  a  trail  of  debris  that  is  now  becoming  a  greater 
concern  for  the  safety  of  future  space  missions.  In  November  of  1995,  the  White  House 
published  the  Interagency  Report  on  Orbital  Debris  [5]  through  the  Office  of  Science 
and  Technology  Policy,  which  is  a  technical  assessment  of  the  orbital  debris  that  has 
accumulated  in  space.  This  report  brought  forth  a  new  sense  of  urgency  to  the  space 
debris  problem  in  the  US. 

A  few  years  later,  the  international  community  looked  forward  to  the  particular 
dangers  of  all  the  debris  left  in  space  and  the  United  Nations  released  a  technical 
report  in  1999  highlighting  some  of  their  concerns  for  the  possible  outcomes.  Af¬ 
ter  many  years  of  conferencing  and  technical  consulting,  the  international  forum  of 
governmental  bodies  called  The  Inter-Agency  Space  Debris  Coordination  Commit¬ 
tee  (IADC)  published  the  IADC  Space  Debris  Mitigation  Guidelines  in  September  of 
2007,  where  they  outline  best  practices  for  the  space  community  to  follow.  These  will 
reduce  the  probability  of  an  undesirable  collision  of  space  debris  with  an  operating 
satellite. 

Although  the  National  Aeronautics  and  Space  Administration  (NASA)  had  been 
looking  at  the  debris  problem  since  at  least  the  late  1980s,  they  finally  made  manda¬ 
tory  some  of  the  widely  accepted  practices  in  the  last  couple  of  years  with  the  publica- 
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tion  of  the  NASA  Technical  Standard  8719.14  (STD  8719.14),  the  NASA  Procedural 
Requirements  8715.6  (NPR),  and  the  NASA  Handbook  8719.14  (NHDBK). 


2.1.1  Science  and  Technology  Council  Update  Report.  During  1995,  the  Na¬ 
tional  Science  and  Technology  Council  (NSTC)  reviewed  the  US  Government’s  1989 
Interagency  Report  on  Orbital  Debris  to  update  its  findings  and  recommendations,  as 
well  as  to  depict  the  progress  and  comprehension  of  the  orbital  debris  environment. 

Later  in  November  they  released  the  1995  Interagency  Report  on  Orbital  Debris, 
where  they  reemphasized  the  growing  problems.  They  noted  that  the  International 
Space  Station  (ISS)  had  taken  steps  to  maximize  protection  from  debris  penetration 
through  implementing  state-of-the-art  shielding,  avoiding  larger  debris,  and  develop¬ 
ing  operational/design  options  that  minimize  the  risk  to  the  station.  [5] 

According  to  the  review,  there  was  a  fairly  clear  picture  of  the  debris  environ¬ 
ment  in  low  earth  orbit.  The  review  also  stated  that  a  particular  concern  was  the 
sustained  rate  of  fragmentation  regardless  of  what  everybody’s  “mitigation  efforts” 
were. 

Another  important  point  was  that  the  NSTC  brought  forth  the  prerequisite  of 
expected  development  of  technical  cooperation  and  consensus  amongst  nations  prior 
to  international  agreement  on  any  regulatory  regimes.  The  concern  for  this  was  the 
possible  competitive  disadvantage  resulting  from  any  unilateral  action  by  the  US  in 
trying  to  diminish  debris  when  nobody  else  is  willing  to  follow  it  or  has  the  technical 
information  to  do  it. 


Table  2.1:  Estimated  debris  population.  [5] 


Size 

Number  of  Objects  %number 

%Mass 

>10  cm 

8,000 

0.02% 

99.93% 

1-10  cm 

110,000* 

0.031% 

0.035%* 

0.1-1  cm 

35,000,000* 

99.67%* 

0.035%* 

Total 

35,117,000* 

100.0%* 

2,000,000  kg# 

*  statistically  estimated  values 


#  calculated  value  from  reported  data 
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Table  2.2:  Catalogued  objects  by  altitude  ranges.  [5] 


Orbit  Type 

LEO 

MEO 

GEO 

Other 

Total 

Cataloged  Objects 

55747 

134 

601 

1447 

7929 

The  estimated  debris  numbers  as  of  November  1995  are  described  by  Table  2.1. 
It  is  important  to  note  that  at  the  time  of  this  report  there  were  a  total  of  about 
8,000  catalogued  objects  in  space  depicted  in  Table  2.2 

2.1.2  United  Nations  Updates  Their  Report.  Years  later,  the  United  Nations 
released  the  Technical  Report  on  Space  Debris  [6]  in  1999.  It  described  the  measuring, 
modeling,  and  mitigation  measures  of  space  debris  for  the  international  community. 

The  Scientific  Technical  Subcommittee,  who  handled  the  report,  articulates  that 
it  was  important  to  have  a  firm  scientific  and  technical  basis  for  future  actions  on  the 
complex  attributes  of  space  debris.  They  later  noted  there  was  already  ongoing  re¬ 
search  in  some  countries  that  helped  the  understanding  of  the  problem  as  a  whole. 
This  research  included  sources  of  debris,  areas  where  the  amount  of  debris  was  in¬ 
creasing,  and  collision  probability  as  well  as  it  effects.  At  last,  this  brought  to  the 
front  the  necessity  to  diminish  the  production  of  debris.  However,  the  report  conceded 
that  there  was  not  even  a  consensus  for  the  definition  of  space  debris.  [6] 

The  report  goes  into  a  summary  of  how  the  community  measures  space  debris. 
An  important  nugget  of  information  was  that  there  was  no  available  information  on 
the  submillimeter  debris  population  above  600  kilometers.  Of  particular  importance 
is  that  there  was  no  information  available  in  the  highest  debris  density  regions  of  800- 
1,000  kilometers  or  the  geostationary  orbit.  When  trying  to  model  the  debris  in  space, 
the  report  used  six  different  debris  environment  models.  Figure  2.1  shows  the  findings 
of  the  different  models  in  spatial  density.  Figure  2.1  has  these  models  showing  that 
the  smaller  the  object  the  more  objects  you  have  in  space  at  a  set  altitude.  However, 
an  even  more  important  graph  is  shown  later  in  the  report  on  Figure  2.2. 


6 


Spatial 

10J 

10-* 


10-5 


10-6 

10-7 

10-* 

10-6 

io-10 

0  200  400  600  800  1  000  1  200  1  400  1  600  1  800  2  000 

Altitude  (km) 

Sources  NASA  (ORDEM96);  DERA  (IDES);  ESA  (MASTER).  CNUCE  (SDM);  and  NAZARENKO 

Figure  2.1:  Model  values  for  current  spatial  density.  [6] 
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Figure  2.2:  Typical  ranges  for  number  of  major  collisions  for  three  scenarios,  1995- 
2095.  [6] 

The  results  of  the  long  term  debris  models  do  not  agree  quantitatively  because  of  differences  in  assumption  and 
initial  conditions.  However,  the  basic  trends  and  tendencies  obtained  by  the  models  agree  qualitatively.  With  the 
continued  implementation  of  debris  mitigation  measures,  the  “business  as  usual”  scenario  could  be  avoided. 


density  (number/knV) 
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This  last  graph  shows  all  models  having  the  same  trends,  an  increase  of  collisions 
in  space.  We  can  safely  deduce  that  debris  will  dramatically  increase  as  collisions 
increase.  The  debris  mitigation  measure  presented  can  be  summarized  as  lessening 
debris  generation  practices  under  normal  operations,  prevention  of  on  orbit  break  ups, 
and  finally  deorbiting  and  reorbiting  of  space  objects. 

According  to  this  report,  the  biggest  mitigation  measure  of  the  space  debris 
problem  was  the  mere  increase  in  awareness  amongst  nations.  The  report  ends  with 
a  somewhat  anticlimactic  tone  when,  in  the  last  paragraph,  it  states  that  “manmade 
space  debris  poses  little  risk  to  . . .  approximately  600  active  spacecraft  now  in  Earth 
Orbit.”  [6] 

2.1.3  IADC  Publishes  Guidelines.  Officially  founded  in  1993,  the  Inter- 
Agency  Space  Debris  Coordination  Committee  used  the  United  Nations  1999  report, 
as  well  as  many  ongoing  proposals  and  space  community  practices,  to  generate  and 
publish  in  September  of  2007  the  IADC  Space  Debris  Mitigation  Guidelines  for  the 
international  community.  The  fundamental  principles  this  publication  follows  are  the 
prevention  of  on-orbit  breakups,  the  removal  of  all  debris  (including  dead  satellites) 
from  useful  regions,  and  the  limiting  of  debris  objects  released  during  normal  opera¬ 
tions.  [8] 

An  interesting  part  of  the  guidelines  are  the  protected  regions  that  must  be 
ensured  for  safe  use.  These  are  the  Low  Earth  Orbit  region  (LEO)  and  the  Geosyn¬ 
chronous  (GEO)  region.  [8]  The  LEO  region,  Region  A  in  Figure  2.3,  consists  of  the 
spherical  region  from  the  surface  of  the  Earth  all  the  way  out  to  2,000  kilometers  in 
altitude. 

The  GEO  region  is  the  a  spherical  shell  described  by  Figure  2.3: 


lower  altitude  =  geostationary  altitude  minus  200  km 
upper  altitude  =  geostationary  altitude  plus  200  km 
-15  degrees  £  latitude  <  +15  degrees 

geostationary  altitude  (Z  geo)  =  35,786  km  (the  altitude  of  the  geostationary  Earth  orbit) 


Figure  2.3:  Protected  regions.  [8] 


2.1.4  NASA  Provides  More  Guidelines.  As  late  as  July  2008,  NASA  ap¬ 
proved  the  Handbook  for  Limiting  Orbital  Debris.  [10]  This  book  is  designed  to  ac¬ 
company  the  previously  published  NASA  STD  8719.14  [9]  and  the  NPR  8715. 6A.  [11] 


Figure  2.4:  Debris  count  as  of  August  2007.  [10] 
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This  latest  update  to  the  guidance  on  NASA’s  orbital  debris  mitigation  is  a 
refinement  and  also  amplifies  some  of  its  applicability  to  resolve  the  problem.  It  goes 
into  detail  and  updates  some  of  the  previous  projections  on  the  amount  of  debris  and 
the  primary  components  depicted  previously  in  Figure  2.4. 

The  regulation  updates  have  mandatory  requirements  on  design  and  implemen¬ 
tation  of  the  mitigation  measures,  hence  showing  a  much  better  commitment  to  help 
with  this  problem.  [5]  It  is  noteworthy  to  see  how  the  debris  of  breakups  in  space 
(whether  from  system  separations,  collisions,  or  unintentional/intentional  detona¬ 
tions,  shown  in  Figure  2.5)  are  really  starting  to  take  charge  of  the  overall  picture: 


anom.  debris 
1.9% 


I  ^  payloads 
25 .2% 


mission -related 
debris 
12.3% 


Figure  2.5:  On-orbit  objects  by  type.  [10] 


As  of  October  2008  there  are  a  total  of  over  12,800  objects  in  space,  with 
less  than  25  percent  of  them  being  payloads,  according  to  Orbital  Debris  Quarterly 
News.  This  currently  amounts  to  over  2,000  tons  of  debris  in  space.  At  the  present 
there  are  models  showing  it  is  no  longer  possible  to  eliminate  the  space  debris  due 
to  the  collision  probability  (that  will  produce  more  debris)  even  if  we  stop  launching 
altogether.  Figure  2.6  clearly  shows  these  projected  increases  in  debris. 
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Figure  1.  Comparison  of  three  different  scenarios.  From  top  to  bottom:  postmission  disposal  (PMD)  only. 
PMD  and  ADR  of  two  objects  per  year,  and  PMD  and  ADR  of  five  objects  per  year,  respectively. 


Figure  2.6:  Debris  projections.  [26] 

2.1.5  Problem  Wrap  Up.  Throughout  the  years,  the  space  community  has 
launched  every  possible  object  into  space,  from  animals  to  sophisticated  satellites, 
all  with  very  important  missions,  but  there  has  not  been  a  very  strong  effort  into 
mitigating  the  space  debris  problem  we  are  now  facing.  It  was  left  alone  to  benefit 
the  advancement  of  technology  and  the  purpose  of  not  adding  another  burden  to  the 
space  acquisition  and  budgets  of  the  world.  Except,  space  operational  regions  are  a 
limited  resource  that  must  be  protected. 

Good  ideas  that  try  to  affect  our  design  and  acquisition  of  new  systems  have 
been  proposed.  It  is  well  understood  that  the  hypervelocities  of  space  objects  make 
them  true  weapons  and  could  eventually  have  disastrous  consequences  to  our  heavily 
technology  dependant  society.  A  better  effort  seems  to  be  brought  forth  by  the  US 
with  the  latest  regulations. 
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2.2  Population  Safety 

The  safety  of  the  population  has  been  put  in  jeopardy  when  reentering  objects 
carrying  dangerous  payloads  don’t  get  incinerated  by  the  extreme  heat  of  reentry  and 
make  it  to  the  ground. 

2.2.1  Skylab  Reentry.  The  very  first  American  Space  Station  was  called 
Skylab  (Figure  2.7).  It  was  launched  on  May  14,  1973.  It  had  many  goals  such  as 
to  prove  humans  can  live  in  space,  to  expand  astronomic  knowledge,  to  detail  earth’s 
resources,  and  to  conduct  experiments,  among  others.  After  the  end  of  the  last  mission 
in  1974  and  after  completing  engineering  tests,  Skylab  was  placed  in  a  stable  orbit 
where  it  would  remain  for  eight  to  ten  years.  Tragically,  the  increase  in  solar  activity 
greatly  increased  the  drag  and  the  station  reentered  back  to  Earth  on  July  11,  1979. 
The  station  might  have  been  over  ninety  tons  in  mass;  therefore,  pieces  of  the  station 
survived  reentry  and,  like  a  display  of  huge  fireworks,  scattered  debris  over  the  Indian 
Ocean  and  Western  Australia.  [35] 


Figure  2.7:  Skylab  workshop,  May  14,  1973.  [35] 
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2.2.2  Cosmos  954  Reentry.  The  dangerous  reentry  on  January  24,  1978  of 
Cosmos  954  occurred  over  Northwest  Canada  and  marked  a  new  sense  of  concern  over 
objects  falling  from  the  darkness  of  space.  It  became  a  high  profile  reentry  for  the  US 
and  Canadian  governments.  The  extreme  importance  was  not  because  the  satellite 
was  designed  to  scan  the  oceans  for  surface  ships,  but  because  it  was  powered  by 
a  small  nuclear  reactor.  The  reactor  was  ninety  percent  enriched  Uranium-235  and 
it  was  “hot  and  alive.”  The  name  “Operation  Morning  Light”  was  the  code  for  the 
joint  operation  the  two  governments  performed  to  clean  up  the  spread  of  radioactive 
fuel.  [36] 

The  path  of  debris  covered  hundreds  of  kilometers.  Dozens  of  radioactive  ob¬ 
jects,  as  well  as  about  four  thousand  particles,  were  spread  over  a  wide  area  ranging 
from  Great  Slave  Lake  to  Saskatchewan  and  Northern  Alberta.  [13]  As  a  settlement, 
Canada  and  the  Union  of  Soviet  Socialist  Republics  agreed  on  a  payment  of  three 
million  Canadian  dollars  for  all  matters  relating  with  the  incident.  [25] 

2.2.3  USA  193  Destruction.  According  to  the  US  Department  of  Defense 
(DoD)  an  experimental  satellite  had  to  be  shot  down  by  a  single  modified  Standard 
Missile-3  (SM-3)  in  February  of  2008.  The  satellite  never  became  operational  after  it 
was  launched  in  December  of  2006.  According  to  the  DoD’s  news  releases,  the  satellite 
could  have  had  over  11,000  kilograms  of  its  mass  survive  the  reentry  which  included 
435  kilograms  of  hazardous  hydrazine  propellant  fuel.  The  successful  destruction  of 
the  satellite  at  10:26pm  EST,  February  20,  2008,  was  widely  publicized.  [1,7,14] 

The  event  brought  a  diplomatic  problem  to  the  forefront.  The  incident  was 
criticized  by  Russia  as  a  test  of  US’s  capability  to  destroy  other  nation’s  satellites 
and  accusing  the  US  of  causing  a  space  arms  race.  It  even  prompted  the  Bush  Ad¬ 
ministration  to  inform  diplomats  that  the  event  should  not  be  seen  under  the  same 
light  as  China’s  destruction  of  its  satellite  the  previous  year.  [2,3,29,31] 
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2.2.4  Other  Famous  Reentries.  At  the  current  time,  there  is  an  average  of 
one  large  piece  of  debris  falling  back  to  the  earth  each  day.  Most  of  them  fall  on  the 
oceans  or  deserts,  but  not  all.  Such  are  the  cases  in  the  following  occurrences  of  a  few 
Delta  II  rocket  bodies.  [24] 


Figure  2.8:  Delta  II  rocket  body  in  Texas.  [24] 


Figure  2.8  . .  is  the  main  propellant  tank  of  the  second  stage  of  a  Delta  2  launch 

vehicle  which  landed  near  Georgetown,  TX,  on  22  January  1997.  This  approximately 
250  kg  tank  is  primarily  a  stainless  steel  structure  and  survived  reentry  relatively 
intact.  ”  [24] 

“On  21  January  2001,  a  Delta  2  third  stage,  known  as  a  PAM-D  (Payload  Assist 
Module  -  Delta),  reentered  the  atmosphere  over  the  Middle  East.  The  titanium  motor 
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Figure  2.9:  Delta  II  Rocket  body  in  Saudi  Arabia.  [24] 

casing  of  the  PAM-D,  weighing  about  10  kg,  landed  in  Saudi  Arabia  about  240  km 
from  the  capital  of  Riyadh.”  (Figure  2.9)  [24] 


Figure  2.10:  Delta  II  Rocket  body  in  South  Africa.  [24] 

“Another  Delta  2  second  stage  reentered  on  21  April  2000  over  South  Africa.  In 
this  incident,  three  objects  were  recovered  along  a  path  nearly  100  km  long:  the  main 
stainless  steel  propellant  tank,  a  titanium  pressured  tank,  and  a  portion  of  the  main 
engine  nozzle  assembly.”  (Figure  2.10)  [24] 
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2.3  Attempt  to  Enforce  a  Safe  Reentry 

2.3.1  Design  Strategy.  Satellite  developers  have  been  under  scrutiny  to  in¬ 
clude  a  better  design  that  will  have  satellites  disintegrate  easier  upon  reentry.  There 
have  been  many  objects  that  have  survived  reentry  because  they  were  made  with 
materials  that  have  high  melting  points.  Some  examples  of  these  materials  are  stain¬ 
less  steel,  titanium,  and  beryllium.  The  satellite  design  must  include  parts  that  have 
lower  melting  points,  and  in  addition,  be  formulated  for  the  system  to  easily  break 
upon  reentry,  making  it  easier  for  the  parts  to  melt  on  their  way  down.  On  the  other 
hand,  if  the  break  up  allows  objects  to  evade  high  melting  points  then  the  design  must 
incorporate  a  method  that  keeps  those  survivable  objects  together  longer  so  that  they 
experience  the  high  melting  temperatures.  Either  one  of  these  answers  must  be  found 
through  the  study  of  breakup  reentry  phenomenon.  [18] 


2.3.2  Tether  Technology.  Studies  have  been  done  to  discover  whether  the 
use  of  electrodynamic  drag  can  be  used  to  decrease  the  orbital  lifetime  of  debris. 
These  studies  have  focused  on  objects  with  an  altitude  of  400  kilometers  and  higher, 
as  shown  in  Table  2.3. 


Table  2.3:  Tether  table  performance.  [17] 


Initial 

Height 

[km] 

Orbit  Inclination 

0° 

25° 

50° 

75° 

D 

E-ORBIT  TI 

HE  [days] 

1400 

170 

220 

325 

EDT 

not 

used 

1300 

140 

185 

280 

1200 

120 

155 

230 

1100 

95 

125 

185 

1000 

70 

95 

140 

375 

900 

55 

70 

110 

280 

800 

45 

55 

80 

200 

700 

30 

40 

55 

140 

600 

20 

30 

40 

80 

500 

15 

20 

25 

40 

400 

10 

15 

15 

20 

Time  to  de-orbit  a  1500  kg  spacecraft  from  a  given  initial  altitude  to 
250  km  with  a  7.5  km  Terminator  Tether™  with  a  mass  of  1%  of  the  S/C. 
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According  to  a  summary  report  of  IADC  AI  19.1,  the  benefits  of  using  tethers 
may  be  significant.  The  first  is  that  spacecraft  may  only  require  between  one  and 
five  percent  of  their  total  mass  on  a  typical  tether  system,  making  them  extremely 
attractive.  A  second  benefit  is  that  tethers  may  decrease  the  de-orbit  time.  A  final 
benefit  is  that  they  reduce  the  Area-Time-Product  (ATP)  of  LEO  orbiting  objects, 
resulting  in  lessening  their  opportunity  for  collision  with  other  objects.  [17] 

Nonetheless,  tethers  are  very  long  and  thin.  Tethers  may  be  between  5  and 
7.5  kilometers  in  length  and  between  0.5  and  2.5  centimeters  wide.  This  increases 
the  possibilities  of  problems  due  to  their  large  increase  in  collision  cross-sectional 
area.  Also,  the  severing  of  the  tethers  can  be  attributed  to  many  culprits  such  as 
meteors,  vibration,  other  debris,  and  even  manufacturing  defects,  with  an  end  result 
of  providing  more  debris.  [17] 

Most  of  the  studies  done  were  under  the  assumption  of  circular  orbits,  where 
tethers  can  be  stable,  but  not  for  any  highly  eccentric  orbits  such  as  Geotransfer 
(GTO).  Furthermore,  the  studies  done  in  IADC  AI  19.1  were  performed  using  different 
particle  flux  models,  creating  the  need  for  a  common  model  where  a  firm  standard 
may  be  defined.  Therefore,  tether  technology  is  in  need  of  further  study  on:  materials 
that  may  survive  impacts,  prolonged  wear  (cuently  are  supposed  to  be  up  there  for 
years),  and  better  space  particle  flux  models  for  collision  probability  studies.  [17] 

2.4  Balloon  Satellites 

The  idea  of  an  object  with  a  balloon  that  is  dozens  of  meters  in  diameter  is 
not  new.  The  satellite  Echo  I  was  launched  in  August  of  1960.  The  mission  of  the 
satellite  was  designed  as  a  passive  communications  deflector.  The  object  weighed  180 
kilograms.  It  was  launched  into  an  orbit  with  a  118  minute  period  and  an  altitude 
of  over  1,500  kilometers.  The  balloon  was  proven  to  work  at  that  altitude  even  with 
micrometeoroids  affecting  it.  [27] 
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Figure  2.11:  The  Echo  I  team  stand  in  front  of  their  balloon.  [21] 

A  second  satellite  balloon  with  an  even  bigger  diameter  of  41  meters  was  launched 
in  January  of  1964  into  an  orbit  with  an  inclination  of  81  degrees.  [4]  One  last  balloon 
was  sent  into  space  in  the  1960s.  The  Pageos  I  was  launched  in  June  of  1966,  but  it 
had  the  diameter  of  Echo  I  and  was  placed  at  an  87  degree  inclination,  with  an  alti¬ 
tude  of  over  4,200  kilometers.  This  satellite’s  mission  was  to  be  a  target  for  geodetic 
purposes.  It  reentered  in  July  of  1975.  [34] 


Figure  2.12:  Pageos  I  inflation  test.  [16] 


18 


2.5  Atmospheric  Unpredictability 

The  atmosphere  is  extreme  in  its  variability;  it  is  in  constant  motion  as  it  is 
affected  by  solar  heating.  As  air  warms  up  it  moves  into  the  higher  parts  of  the  atmo¬ 
sphere.  At  the  same  time,  cooler  air  drops  in  altitude.  This  boiling-type  effect  hap¬ 
pens  vertically  as  well  as  horizontally.  All  this  merging  of  different  air  temperatures 
and  the  rotation  of  the  earth,  coupled  with  the  filtering  of  the  different  atmospheric 
component  weights  (as  altitude  increases  over  120  kilometers)  makes  for  very  erratic 
atmospheric  profiles.  [30, 33] 


Figure  2.13:  Atmosphere  picture  taken  from  the  ISS.  [12] 


There  are  many  different  atmospheric  models  that  may  be  used  when  solving 
for  an  object’s  path  in  earth’s  orbit.  These  models  state  the  atmosphere’s  density  at 
a  certain  location.  The  commonly  used  Standard  1976  atmospheric  density  model  is 
an  updated  version  that  was  originally  published  in  1958  with  updates  in  1962,  1966, 
and  finally  1976.  It  consists  of  tables  that  state  the  temperature,  pressure,  and  mass 
densities  of  the  different  atmospheric  components  as  a  function  of  altitude.  [15,30] 
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Most  atmospheric  density  models  go  through  updates  because  of  the  compli¬ 
cated  wavering  nature  of  the  atmosphere  and  our  continuous  renewing  of  atmospheric 
education.  A  model  that  went  through  the  update  phases  (and  is  still  going  through 
them)  is  the  Mass  Spectrometer  and  Incoherent  Scatter  Radar  (MSIS)  model.  For  this 
reason,  the  family  of  models  for  this  type  came  to  be  known  as  the  MSIS-77,-83,-86, 
with  the  last  one  being  combined  with  the  CIRA-1986.  The  CIRA  models  were  first 
produced  by  the  Committee  on  Space  Research  (COSPAR).  The  latest  MSIS-E  1990, 
where  the  “E”  was  added  because  it  covered  from  the  ground  to  1,000  kilometers  in 
altitude,  is  an  update  that  is  still  very  widely  used.  [22,28]  One  last  permutation  of 
the  MSIS  model  is  the  Navy  Research  Lab  MSISE  2000  (NRLMSISE-00).  This  latest 
model  has  updated  databases  and  uses  inputs  as  solar  flux  and  geomagnetic  heating 
to  report  the  atmospheric  densities.  [23] 

These  never  ending  studies  and  updates  of  atmospheric  density  profiles  clearly 
depict  how  complicated  and  difficult  is  to  even  have  an  accurate  model  of  atmospheric 
density.  As  stated  before,  all  fluctuations  are  directly  related  to  the  sun’s  activities 
which  create  differential  solar  winds  affecting  our  atmosphere.  See  Figure  2.14.  [15,30] 


Solar  wind ,  as  depicted  in  this  artist’s  illustration,  travels  from  the 
Sun  and  envelopes  the  Earth’s  magnetic  field.  High-energy  pulses 
of  solar  wind  from  sunspot  activity  (“solar  bursts”  or  “plasma 
bubbles”)  travel  from  the  Sun  to  the  Earth  at  speeds  exceeding  500 
miles  per  second.  The  pulses  distort  the  Earth’s  magnetic  field  and 
produce  ge magnetic  storms  that  disrupt  the  Earth’s  environment. 


Figure  2.14:  Solar  wind  disrupting  earth’s  magnetic  field.  [30] 
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The  sun  has  its  own  disturbances  and  the  activities  seem  to  have  a  cycle  of 
about  eleven  years.  When  the  sun  is  most  active,  it  is  said  to  be  in  “solar  max”  and 
when  the  activity  slows  down  it  is  called  “solar  min.”  These  changes  in  activity  have 
marked  effects  on  our  atmosphere.  During  these  activities  the  sun  releases  tremendous 
amounts  of  energy  and  mass,  producing  strong  winds.  These  winds  change  the  earth’s 
atmosphere  making  it  compressed  and  changing  the  drag  satellite  systems  perceive. 
Thus,  the  sun’s  activities  have  a  great  effect  on  the  true  atmospheric  drag  that  objects 
orbiting  the  earth  “feel”  and  our  current  models  can  only  try  to  predict  what  an 
orbiting  path  will  look  like.  [30] 


Figure  2.15:  Solar’s  unusual  activities  during  Solar  Min.  [32] 
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2.6  Square  Cube  Law 

The  square  cube  law  was  first  introduced  by  Galileo  Galilei  in  his  writings  called 
Two  New  Sciences.  He  discusses  how  an  object  that  increases  proportionally  in  size 
will  not  have  the  same  proportional  increase  in  strength.  This  phenomenon  has  to 
do  with  the  area  to  volume  ratio,  which  is  commonly  used  as  the  area  to  mass  ratio, 
since  it  is  assumed  to  keep  a  constant  density.  Therefore,  when  you  increase  the  size 
of  an  object  the  surface  area  is  increased  by  the  square  of  its  length  wile  the  mass  is 
increased  by  the  cube  of  its  length.  This  means  the  object  will  become  heavier  much 
faster  than  the  surface  area  it  will  cover.  [19, 20] 

This  can  be  seen  by  the  simple  illustration  of  a  cube.  A  cube  with  a  length  L 
will  have  the  surface  area  of  SA  —  6  *  L2  and  a  volume  of  V  —  L3.  Therefore,  if 
the  length  of  the  object  is  increased  by  an  order  of  magnitude,  the  cube  will  have  a 
length  of  L  =  10,  a  surface  area  equal  to  SA  —  6  *  (10  *  10)  =  600,  and  a  volume  of 
V  =  (10  *  10  *  10)  =  1000.  This  shows  that  the  object  increased  its  length  by  a  factor 
of  ten,  its  surface  area  by  a  factor  of  one  hundred,  and  its  volume  (or  mass)  by  one 
thousand.  [19,20] 


10  L 


Cube  A 


Cube  B 


Figure  2.16:  Increase  in  length  demonstration.  [19] 
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This  law  explains  many  phenomenon,  including  the  size  of  animals,  how  they 
cope  with  the  need  to  keep  warm,  and  even  their  ability  to  move  around.  An  elephant 
is  not  capable  of  jumping  because  its  weight  prohibits  it,  but  an  ant  falling  from  a 
skyscraper  will  survive  and  go  on  to  carry  fifty  times  its  own  weight.  Even  biological 
cells  stay  at  the  microscopic  size  levels,  because  the  high  efficiency  required  to  transfer 
nutrients  and  waste  across  their  membranes  occurs  at  those  sizes.  Consequently,  cells 
just  multiply,  keeping  their  high  surface  area  to  mass  ratio  instead  of  growing.  [19,20] 


Notice  that  the  shape  of  the  mouse  is  as  close  as  it  can  be  to  that  of  a  sphere,  a  shape  that  produces  the  minimum  amount 
of  surface  area.  Now  notice  that  the  elephant's  skin  is  covered  with  wrinkles  thus  increasing  the  elephant's  surface  area. 
Also  notice  that  the  mouse  is  like  nearly  all  mammals  in  being  covered  with  a  thick  fur  coat  while  the  elephant  is  nearly 
void  of  all  hair.  The  small  mouse  is  trying  to  maintain  heat  while  the  large  elephant  is  trying  to  lose  it. 

Figure  2.17:  Drastic  surface  area  comparison.  [19] 

This  thesis  will  use  the  fact  that  when  inflating  a  balloon,  the  surface  area 
increases  by  the  square  of  the  radius  while  at  the  same  time  decreasing  the  actual 
density  of  the  object  as  a  whole.  Although  a  balloon  is  the  geometric  figure  with  the 
smallest  surface  area  to  volume  ratio,  it  is  also  the  most  consistent  when  viewed  from 
any  angle.  This  means  that  it  does  not  have  to  be  oriented  for  it  to  use  its  maximum 
cross-sectional  area  when  it  is  being  used  as  the  agent  to  increase  drag. 
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III.  Methodology 


All  simulations  for  this  thesis  were  performed  using  the  Satellite  Tool  Kit  (STK) 
software  produced  by  Analytical  Graphics,  Inc  (AGI).  This  entire  project  will 
use  the  International  System  of  Units  (SI). 

STK  is  a  commercial  off  the  shelf  analysis  software  that  may  be  used  to  depict 
scenarios  in  a  very  clean  and  clear  format.  The  software  is  capable  of  using  three- 
dimensional  visualization  and  using  many  kinds  of  integration  algorithms  to  answer 
problems.  It  solves  fundamental  location  and  inter-visibility  problems  associated  with 
land,  sea,  air,  and  space  scenarios.  The  overall  purpose  for  each  of  the  procedures  will 
be  stated  first,  followed  by  the  specific  simulation  set  up  including  parameters  used 
along  with  the  rationality  behind  them. 

3.1  Initial  Set  Up  and  First  Study:  Variability  in  Reentry  Prediction 
Due  to  the  Atmosphere 

The  atmosphere’s  unpredictability  is  an  important  factor  that  must  be  taken 
into  account  in  any  kind  of  analysis  dealing  with  an  object’s  reentry.  To  explore  the 
variability  in  reentry  prediction  due  to  atmosphere  variation,  a  particular  orbit  was 
propagated  using  different  atmospheric  models.  Furthermore,  two  additional  orbits 
were  used  in  the  simulations  to  account  for  the  reentry  variability  in  different  orbit 
types. 


3.1.1  Propagator.  All  orbits  were  propagated  using  the  High  Precision 
Orbit  Propagator  in  STK.  The  program  has  the  option  of  setting  an  atmospheric 
model  from  ten  different  available  models  that  have  been  used  through  the  years. 
Thus,  in  order  to  get  the  best  picture  of  what  the  reentry  could  look  like,  all  ten 
atmospheric  models  were  used.  The  mere  fact  that  there  are  ten  models  to  choose 
from  shows  the  unpredictability  of  the  atmosphere. 

The  available  atmospheric  density  models  in  STK  are  the  following: 

1.  1976  Standard 
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2.  CIRA  1972 


3.  Harris-Priester 

4.  Jacchia  1960 

5.  Jacchia  1970 

6.  Jacchia  1971 

7.  Jacchia-Roberts 

8.  MSIS  1986 

9.  MSISE  1990 

10.  NRLMSISE  2000 

For  the  rest  of  the  study,  especially  throughout  the  plots  and  analysis,  they  will 
mostly  be  referred  to  by  the  number  they  were  assigned  in  this  section  instead  of  their 
whole  name  in  order  to  add  clarity  to  the  visual  results. 

3.1.2  Reentry  Criterion.  A  reentry  criterion  must  be  set  for  the  study.  The 
reason  for  this  is  because  several  atmospheric  models  are  not  capable  of  modeling 
below  90  kilometers.  Air  Force  Space  Command  considers  an  object  reentered  at  120 
kilometers.  Using  these  two  criterion,  100  kilometers  was  used  as  the  point  where  an 
object  will  be  considered  reentered.  Even  though  the  reentry  criterion  was  set,  all 
simulations  were  executed  all  the  way  to  the  lowest  atmospheric  altitude  limit  set  by 
their  correspondent  atmospheric  model. 

3.1.3  Simulated  Orbits.  Three  different  orbital  element  sets  were  chosen 
for  the  study  in  order  to  ensure  a  wide  enough  number  of  factors  are  covered.  The 
particular  element  sets  can  be  found  in  Appendix  A. 

1.  The  first  propagated  orbital  element  set  was  from  a  piece  of  debris  that  was  left 
from  the  destruction  of  USA  193.  This  orbit  has  a  58.5  degree  inclination  with 
an  eccentricity  of  .01687  and  a  period  of  91.4  minutes.  The  specific  numbers 
used  to  set  up  STK  are  on  Figure  B.l  in  Appendix  B. 
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2.  The  second  orbital  element  set  used  was  from  the  ISS  and  has  a  51.7  degree 
inclination  with  an  eccentricity  of  .00100  and  a  period  of  91.6  minutes.  The 
specific  numbers  used  to  set  up  STK  are  in  Figure  B.2  in  Appendix  B. 

3.  The  third  element  set  used  was  from  a  satellite  named  Terrasar-X  in  a  polar 
orbit.  It  has  a  97.4  degree  inclination  with  an  eccentricity  of  0.00014  and  a  period 
of  94.8  minutes.  The  specific  numbers  used  to  set  up  STK  are  in  Figure  B.3  in 
Appendix  B. 


3.1.4  Object  Parameters.  To  get  uniform  results,  the  same  satellite  dimen¬ 
sions  were  used  for  all  simulations.  The  object  has  the  default  mass  of  1,000  kilograms, 
with  an  area  to  mass  ratio  of  0.1  and  an  inertial  matrix  that  looks  like  Table  3.1. 


Table  3.1:  Satellite’s  inertial  matrix. 


X 

Y 

Z 

X 

4500  kg  m2 

0 

0 

Y 

0 

4500  kg  nr2 

0 

Z 

0 

0 

4500  kg  m2 

3.1.5  Trajectory  Data.  Once  the  orbit  was  propagated,  STK  provided  po¬ 
sitional  data  in  Latitude,  Longitude,  Altitude  (LLA)  tables  for  the  entire  satellite 
trajectory.  An  LLA  table  includes  the  previously  mentioned  position  parameters  as 
well  as  corresponding  times  and  rates  of  change.  STK  is  able  to  provide  these  tables  in 
a  spreadsheet  format.  Using  Microsoft  Excel,  all  reentry  path  altitudes  of  the  object 
were  plotted  against  time  to  explore  their  reentry  variability. 

3.1.6  Chosen  Atmospheric  Model.  Atmospheric  model  number  9,  which  is 
the  MSISE  1990,  was  chosen  to  be  used  as  the  effective  reference  reentry  path  for  the 
rest  of  the  simulations.  I11  other  words,  the  reference  reentry  path  will  be  used  as  if 
it  was  the  actual  path  the  satellite  would  fallow  to  reentry.  The  rationalization  for 
this  was  that  the  model  predictions’  were  the  closest  to  the  median  out  of  all  models 
predictions. 
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In  addition,  this  atmospheric  model  is  capable  of  propagating  to  zero  altitude. 
The  STK  program  has  an  option  to  increase  the  drag  by  an  input  scale  at  an  input 
time.  This  option  was  used  for  the  following  set  of  simulations.  The  goal  for  this  part 
of  the  experiment  is  to  investigate  how  a  drag  increase  by  a  particular  factor  affects 
the  overall  reentry  time. 

3.2  Second  Study.  Drag  Scale  Increase  Study 

The  drag  in  the  propagator  was  increased  by  a  factor  of  10.  The  drag  increment 
was  set  for  4  separate  position  points  through  the  orbit.  Once  the  drag  was  increased, 
the  orbit  was  propagated  all  the  way  to  whatever  atmospheric  limitation  was  set  by 
the  model.  For  every  single  simulation,  STK  provided  data  in  spreadsheets  that  were 
later  used  to  plot  the  reentry  path  altitudes  as  a  function  of  time. 

1.  The  first  set  of  simulations  for  this  study  were  using  the  element  set  of  the 
debris  from  USA  193.  The  drag  scale  increase  was  set  at  the  beginning  of  the 
day  of  the  reference  reentry,  regardless  of  the  time  left  to  reentry.  To  observe 
how  much  the  total  time  to  reentry  would  change,  the  next  3  drag  increments 
were  set  24  hours  apart. 

One  more  set  variable  was  the  point  in  the  orbit  where  the  drag  increments 
happen.  In  order  to  detect  any  blatant  effect  due  to  the  location  of  the  drag 
increment,  the  drag  increments  were  alternated  between  perigee  and  apogee. 
The  first  drag  increase  happened  at  perigee  which  was  closer  to  the  reference 
reentry  time. 

2.  The  second  set  of  simulations  was  on  the  element  set  from  the  ISS.  This  first 
drag  scale  increase  also  happened  at  the  beginning  of  the  day  of  reentry.  Again, 
to  observe  how  the  reentry  time  between  the  drag  increase  and  final  reentry 
would  change,  the  next  drag  increments  were  set  even  further  apart,  at  4  days. 

3.  The  last  set  of  simulation  was  on  the  Terrasar-X.  The  first  drag  scale  increase 
was  also  at  the  beginning  of  the  day  of  reentry.  This  time,  in  order  to  observe 
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the  change  in  time  left  to  reentry,  the  rest  of  the  drag  increments  were  about 

24  hours  apart. 

3.2.1  Greater  Scale  Increase.  To  briefly  explore  what  a  larger  increase  in 
the  drag  scale  would  do  to  the  time  of  reentry,  the  scale  was  increased  to  100  and  the 
same  procedures  as  above  were  followed. 

3.3  Third  Study.  Single  A/M  Increase  at  Several  Points  in  Orbital 

Path 

This  thesis  proposes  to  increase  the  drag  force  exerted  on  the  object  through 
the  amplification  of  its  surface  area.  It  is  a  more  intuitive  way  of  thinking  about  the 
reentry  process.  The  greater  the  surface  area,  the  greater  the  acceleration  (a^)  due 
to  drag  the  satellite  will  experience  according  to  the  basic  equation: 

ad  =  -\cdp—v2  (3.1) 

2  m 

[37, 38]  where  the  symbols  represent  the  following: 

•  Cd  —  Object’s  Coefficient  of  Drag 

•  p  —  Atmospheric  Density 

•  A  —  Object’s  Cross-sectional  Area  Perpendicular  to  Velocity 

•  m  —  Object’s  Mass 

•  v  =  Object’s  Velocity 

Thus,  when  computing  drag  change,  it  is  easier  to  visualize  changing  a  specific 
parameter  in  the  simulation.  This  factor  would  be  the  area  to  mass  ratio  (A/M)  and 
STK  has  the  capability  to  adjust  it. 

The  area  used  in  the  equation  is  the  cross  sectional  area  perpendicular  to  the 
velocity  vector  of  the  orbiting  object.  Consequently,  for  the  rest  of  the  simulations, 
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the  capability  of  modifying  the  area  to  mass  ratio  will  be  used  instead  of  adjusting 
the  drag  by  a  certain  factor. 

This  next  set  of  simulations  will  look  directly  into  the  variability  of  the  reentry 
time  due  to  the  atmospheric  models.  The  A/M  ratio  will  be  increased  from  the  default 
of  0.1  to  1.0  for  all  the  different  atmospheric  models.  The  results  will  be  compared 
against  the  reference  orbit  solved  by  the  MSISE  1990  with  a  0.1  A/M  ratio. 

The  capability  of  STK  to  increase  the  drag  by  a  certain  factor  had  the  option 
of  inputting  the  specific  time  when  the  drag  would  increase.  On  the  other  hand,  the 
capability  to  manipulate  the  A/M  ratio  is  set  for  the  whole  propagation  and  cannot 
be  made  to  change  at  a  certain  time.  This  prompted  the  use  of  an  element  set  (that 
was  solved  for)  from  the  reference  orbit.  Therefore,  at  the  chosen  time  for  the  increase 
in  A/M  ratio,  the  element  set  from  the  reference  orbit  was  used  as  the  epoch  of  the 
orbit  and  the  simulations  began  from  that  point  to  final  decay. 

For  every  single  simulation,  STK  provided  data  in  spreadsheets  that  were  later 
used  to  plot  the  reentry  path  altitude  against  time. 

1.  The  USA  193  debris  orbital  element  was  used  first.  The  first  increase  in  A/M 
ratio  was  a  little  more  than  2  days  and  9  hours  before  the  reference  reentry  time. 
The  second  A/M  ratio  increase  was  done  24  hours  prior  to  reference  reentry. 
The  last  ratio  increase  was  set  10  hours  prior  to  reference  reentry. 

2.  The  next  simulation  was  the  ISS  orbital  element  set.  The  first  increase  in  A/M 
ratio  was  7.5  days  before  the  reference  reentry  time.  The  second  A/M  ratio 
increase  was  done  24  hours  prior  to  reference  reentry.  The  last  ratio  increase 
was  set  10  hours  prior  to  reference  reentry. 

3.  The  last  simulation  was  the  Terrasar-X  element  set.  The  first  set  to  be  propa¬ 
gated  was  3  days  before  reentry.  The  second  A/M  ratio  increase  was  24  hours 
prior  to  reentry  and  the  last  ratio  increase  was  10  hours  before  reentry. 
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3.4  Fourth  Study  :  A/M  Manipulation  at  a  Single  Point  in  Orbital  Path 

The  last  set  of  simulations  was  focused  on  briefly  exploring  how  the  A/M  ratio 
affected  the  reentry  prediction  time  as  it  was  manipulated.  The  only  element  set  to 
be  propagated  for  this  part  of  the  study  is  the  original  USA  193  debris  orbital  element 
set.  The  atmospheric  model  used  throughout  will  be  the  MSISE  1990.  Only  the  A/M 
ratio  will  be  changed.  Eleven  different  ratios  were  chosen,  spread  between  0.001  and 
1.000.  The  reentry  altitudes  were  plotted  on  3  different  plots  to  study  their  reentry 
behavior  more  clearly. 
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IV.  Simulation  Results  &  Analysis 


The  four  different  studies  show  firm  results  of  the  effects  of  atmospheric  unpre¬ 
dictability,  drag  increase,  and  the  amplification  of  the  object’s  surface  area. 

4-1  Atmospheric  Unpredictability 

Here  are  the  simulations  results  and  discussions  of  the  study  on  the  variability 
of  reentry  predictions  due  to  atmospheric  instability. 

4-1.1  Debris  Reentry  Prediction.  The  propagation  path  of  the  USA  193 
debris  element  set  using  all  atmospheric  models  is  shown  in  Figure  4.1: 


Time  (days) 


Figure  4.1:  Debris  orbit  reentry  predictions  by  ten  atmospheric  models. 
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The  atmosphere’s  variability  showed  it  has  a  big  impact  on  the  reentry  time 
prediction.  The  period  between  the  epoch  of  the  element  set  until  the  reference 
reentry  time  was  almost  116  hours.  This  means  that  the  element  set  was  propagated 
almost  5  days  before  reentry.  The  result  was  a  reentry  prediction  spread  of  over  47 
hours  among  the  ten  different  atmospheric  models. 

Taking  out  the  predictions  that  are  not  with  in  a  few  hours  of  the  reference 
prediction  (3  outliers)  from  the  group,  the  time  spread  dropped  down  to  almost  5.5 
hours.  Still,  a  323  minute  spread  amounts  to  over  3.5  revolutions  around  the  earth, 
making  the  predictions  useless.  The  results  are  summarized  in  Figure  4.2: 
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Figure  4.2:  Summarized  predicted  reentry  times  for  Debris  orbit. 
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4-1-2  ISS  Orbit  Reentry  Prediction.  The  propagation  of  the  ISS  orbit  in  all 
ten  atmospheric  models  provided  the  reentries  show  in  Figure  4.3: 


Figure  4.3:  ISS  orbit  reentry  predictions  by  ten  atmospheric  models. 

The  time  from  initial  orbit  propagation  until  reference  reentry  was  approxi¬ 
mately  13  days  or  almost  316  hours.  The  actual  time  span  of  all  reentry  predictions 
was  over  6  days,  or  more  than  148  hours,  shown  by  the  summarized  results  Figure  4.4: 
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Figure  4.4:  Summarized  predicted  reentry  times  for  ISS  Orbit. 

Again,  for  this  simulation,  there  seem  to  be  3  predictions  that  are  outliers. 
These  are  the  3  earliest  ones,  which  are  separated  by  more  than  a  few  hours  from  the 
reference  reentry.  When  taking  these  points  out  of  the  spread  the  window  is  reduced 
to  almost  16  hours.  This  equals  to  over  10  revolutions  around  the  earth,  making  the 
predictions  hollow. 

4-1.3  Terrasar-X  Prediction.  The  last  element  set  propagated  for  the  study 
was  that  of  the  Terrasar-X  satellite.  There  were  some  extreme  results  for  these  sim¬ 
ulations  as  depicted  in  Figure  4.5  and  Figure  4.6: 
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Figure  4.5:  Terrasar-X  orbit  reentry  predictions  by  seven  atmospheric  models. 

In  view  of  the  fact  that  this  is  a  very  different  orbit  from  those  previously  used, 
the  reference  prediction  is  more  than  6  months  apart  from  its  initial  propagation 
point.  The  reason  for  the  extreme  time  to  reenter  is  the  orbits’  low  eccentricity  and 
the  much  higher  altitude  of  over  500  kilometers. 

The  long  orbit  propagation  only  increased  the  effect  of  the  atmospheric  vari¬ 
ability  in  the  reentries,  increasing  the  prediction  spread  to  almost  79  days.  To  see  the 
reentries  more  clearly,  the  graph  was  split  amongst  the  2  groups. 
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Figure  4.6:  Terrasar-X  orbit  reentry  predictions  by  three  atmospheric  models. 

The  first  group  of  predictions  (shown  in  Figure  4.5)  are  the  predictions  that  are 
closer  to  each  other.  The  second  group  of  predictions  in  Figure  4.6  were  considered  to 
be  outliers  and  are  the  ones  that  contribute  the  most  to  the  prediction  time  spread. 

When  only  considering  the  first  group  of  predictions  for  the  analysis,  the  reentry 
spread  dropped  down  to  over  3.5  days.  Again,  this  84  hour  prediction  spread  is  of 
almost  no  particular  use.  The  specific  numbers  are  summarized  in  the  Figure  4.7: 
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Figure  4.7:  Summarized  predicted  reentry  times  for  Terrasar-X  orbit. 


4-1-4  Atmospheric  Unpredictability  Conclusion.  These  simulations  have 
confirmed  the  impracticality  of  any  reentry  prediction  that  is  more  than  a  few  days 
into  the  future  due  to  the  unpredictability  of  the  atmosphere.  As  a  result,  the  study 
showed  a  variability  increase  in  reentry  prediction  as  the  actual  reentry  was  further 
into  the  future. 
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\.2  Drag  Increase  by  a  Factor 

The  STK  software  program  provides  the  capability  to  increase  drag  by  a  certain 
factor  at  a  certain  point  in  time.  This  feature  was  used  to  explore  the  effects  of 
increased  drag,  as  in  the  case  of  solar  panel  deployments,  at  different  points  in  the 
orbital  path.  This  study  focuses  on  the  time  reentry  prediction  only  and  not  the 
location  reentry  prediction. 

4-2. 1  Debris  Orbit  Reentry.  Figure  4.8  shows  the  dramatic  change  the  object 
experiences  in  its  reentry  path  when  the  drag  is  scaled  up  by  a  factor  of  ten.  The 
drag  is  increased  at  a  single  point  and  propagated  until  reentry.  This  is  done  for  four 
different  points  (not  on  the  same  simulation).  Figure  4.8  shows  all  five  simulations: 


Figure  4.8:  Debris  reentry  predictions  due  to  drag  increase  at  several  points  in  the 
orbit  path. 


38 


Each  drag  increase  is  set  at  the  beginning  of  the  calendar  day.  In  the  first  drag 
scale,  the  object  accelerated  and  took  just  over  9  hours  to  reenter  with  the  reference 
reentry  prediction  almost  4  days  into  the  future.  Thus,  the  same  thing  happens  for 
the  rest  of  the  points.  It  is  important  to  note  that  the  graph  shows  an  increase  in  the 
objects’  acceleration  the  closer  to  reentry  the  drag  increase  is  set.  This  is  depicted  by 
the  top  graph  in  Figure  4.9: 
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Figure  4.9:  Predicted  reentry  times  for  Debris  orbit;  Original  vs.  New  Time  to  Decay 
plot;  New  over  Original  Time  to  Decay  Ratio  plot. 
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In  fact,  the  ratio  of  the  new  time  to  reentry  over  the  reference  time  to  reentry 
increases  as  you  get  closer  to  reentry.  This  means  the  closer  the  object  is  to  reentry 
the  more  the  new  time  to  reentry  approaches  the  original  time  to  reentry.  This  is 
expected  since  it  should  approach  one  at  the  point  of  reentry.  Still  the  ratio  stays  in 
the  teens  for  the  4  chosen  points. 

To  see  how  the  reentries  behaved  with  an  extreme  drag  increase,  the  factor 
by  which  the  drag  was  increased  was  set  to  100  and  path  changes  are  depicted  in 
Figure  4.10. 


Time  (days) 


Figure  4.10:  Debris  reentry  predictions  due  to  extreme  drag  increase  at  several  points 
in  the  orbit  path. 


40 


It  is  obvious  that  with  a  drag  increase  of  this  magnitude  the  object  was  made 
to  come  in  within  1  revolution. 


Figure  4.11  shows  the  new  time  predictions  to  reentry  and  how  the  ratio  between 
the  new  prediction  and  the  old  prediction  did  not  increase  as  dramatically  because  of 
the  already  small  ratio  values: 
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Figure  4.11:  (Second)  Predicted  reentry  times  for  Debris  orbit;  Original  vs.  New 
Time  to  Decay  plot;  New  over  Original  Time  to  Decay  Ratio  plot. 
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4-2.2  ISS  Orbit  Reentry.  This  second  set  of  simulations  focuses  on  exploring 
the  change  in  the  time  left  to  reentry  by  increasing  the  time  in  between  points  of  drag 
increases  from  24  hours  in  the  previous  set  of  simulations  to  96  hours.  Figure  4.12 
illustrates  the  reentry  paths: 


Figure  4.12:  ISS  reentry  predictions  due  to  drag  increase  at  several  points  in  the  orbit 
path. 

The  first  drag  scale  increase  was  set  over  12.5  days  prior  to  reentry  and  the  time 
between  drag  increase  points  was  set  to  4  days.  Figure  4.13  shows  the  same  trends 
from  the  previous  simulations: 
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Figure  4.13:  Predicted  reentry  times  for  ISS  orbit;  Original  vs.  New  Time  to  Decay 
plot;  New  over  Original  Time  to  Decay  Ratio  plot. 

The  interesting  change  was  that  the  ratio  of  the  new  time  to  reentry  over  the 
old  time  to  reentry  showed  a  little  stagnation  until  the  third  point.  The  ratios  are  still 
in  the  10  to  11  percent  range  of  the  original  time  left  to  reentry.  Figure  4.14  shows 
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the  reentry  paths  resulting  from  an  increase  in  drag  scale  of  100: 
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Figure  4.14:  ISS  reentry  predictions  due  to  extreme  drag  increase  at  several  points  in 
the  orbit  path. 


Once  more,  the  extreme  increase  in  drag  pulled  the  reentries  to  within  4  hours 
from  the  drag  increase.  The  following  summary  shows  the  last  one  making  it  reenter 
in  36  minutes: 
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Figure  4.15:  (Second)  Predicted  reentry  times  for  ISS  orbit;  Original  vs.  New  Time 
to  Decay  plot;  New  over  Original  Time  to  Decay  Ratio  plot. 


Another  observation  is  that  in  the  extreme  scale  increase  of  100,  the  ratios  of 
time  left  to  reentry  were  in  the  single  percent  digits,  (see  Figure  4.15)  and  in  the  scale 
increase  of  10,  shown  in  Figure  4.13,  the  time  left  to  reentry  ratios  were  in  the  low 
teens. 
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4-2.3  Terras ar-X  Orbit  Reentry.  The  last  set  of  simulations  for  this  study 
explore  the  effects  of  the  same  drag  scale  increases  used  in  the  previous  2  sets.  This 
type  of  orbit  will  show  what  kind  of  an  effect  the  increase  in  inclination  will  have  on 
the  patterns  shown  before  by  the  increases  in  drag. 

For  the  Terrasar-X  element  set,  the  first  scale  increase  was  set  3.5  days  prior  to 
reference  reentry  and  the  rest  were  spread  out  in  24  hour  intervals. 


Figure  4.16:  Terrasar-X  reentry  predictions  due  to  drag  increase  at  several  points  in 
the  orbit  path. 


Figure  4.16  shows  the  same  behavior  as  the  2  previous  sets  of  simulations,  where 
new  time  to  reentry  drops  as  you  get  closer  to  reference  reentry  from  almost  9  hours 
to  96  minutes,  depicted  in  Figure  4.17: 
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Figure  4.17:  Predicted  reentry  times  for  ISS  orbit;  Original  vs.  New  Time  to  Decay 
plot;  New  over  Original  Time  to  Decay  Ratio  plot. 


One  last  time  the  drag  scale  was  increased  to  100  and  the  same  conduct  from 
the  previous  simulations  was  observed.  The  reentry  path  dramatically  accelerated 
towards  reentry. 
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Figure  4.18:  Terrasar-X  Reentry  predictions  due  to  extreme  drag  increase  at  several 
points  in  the  orbit  path. 


Similar  results  were  observed  as  depicted  in  Figure  4.19  and  again  the  time  to 
reentry  ratios  in  the  bottom  graph  of  Figure  4.19  show  them  to  be  in  the  single  digit 
percentile. 
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Figure  4.19:  (Second)  Predicted  reentry  times  for  ISS  orbit;  Original  vs.  New  Time 
to  Decay  plot;  New  over  Original  Time  to  Decay  Ratio  plot. 


4-2-4  Drag  Increase  Conclusion.  This  study  has  established  making  an 
object  reenter  through  the  use  of  drag  is  achievable  with  the  appropriate  increase  in 
drag.  It  was  observed  that  the  ratio  of  the  new  time  to  reentry  over  the  old  time  to 
reentry  seemed  to  be  inversely  proportional  to  the  scale  of  drag  increase.  Also,  the 
drag  increase  shows  a  linear  relationship  to  the  time  left  to  reentry. 
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4-3  Reentry  Prediction  with  a  Single  Increase  to  A/M  at  Several  Points 
in  Orbit  Path 

In  order  to  investigate  more  closely  what  the  increase  in  surface  area  would  do 
to  the  reentry  path  of  an  object  the  A/M  ratio  was  manipulated  in  STK.  All  3  orbits 
previously  used  in  the  studies  were  used  for  this  study.  There  were  3  points  chosen 
in  the  path  of  the  reference  reentry  where  the  A/M  ratio  was  increased  from  0.1  to 
1.0.  The  furthest  point  away  from  the  reentry  time  was  simulated  first  and  all  ten 
atmospheric  models  were  used  for  the  plotted  propagations. 

4-3.1  Debris  Variability  Change.  The  first  element  set  propagated  was  the 
Debris  set.  The  first  set  of  simulations  increased  the  A/M  ratio  at  the  57  hour  point 
before  the  reference  reentry.  Figure  4.20  shows  the  reentry  prediction  paths,  with  the 
larger  A/M  ratio,  by  all  atmospheric  models  against  the  reference  reentry.  It  displays 
a  rapid  descent  in  altitude  with  all  the  reentry  predictions  close  to  each  other. 

4-3. 1.1  First  A/M  Increase.  In  Figure  4.21,  the  new  reentry  predic¬ 
tions  are  shown  without  the  reference  reentry  to  hone  in  on  their  variability  spread. 
The  time  spread  between  the  earliest  prediction  and  the  latest  one  was  just  over  2 
hours.  It  is  starting  to  become  clear  that  the  time  prediction’s  uncertainty  is  decreas¬ 
ing  greatly.  The  original  57  hours  left  to  reentry  decreased  to  about  6  hours  and  the 
original  reentry  time  prediction  spread  of  2  days  (from  Figure  4.2),  has  dropped  to  a 
forced  reentry  spread  of  2  hours. 

Furthermore,  the  fact  that  the  reentry  time  has  dropped  to  nearly  2  hours  indi¬ 
cates  the  reentry  locations  are  a  bit  more  than  a  revolution  apart.  Thus,  Figure  4.22 
shows  4  reentry  points  off  the  trend-line  but  6  of  them  start  to  expose  a  gain  in  pre¬ 
cision.  At  this  point  you  are  able  to  predict  the  location  of  the  reentry  within  a  ring 
around  the  earth  where  you  essentially  would  want  to  land  it  or  not. 
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Figure  4.20:  Debris  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  57  hours  prior  to  original 
reentry. 
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Figure  4.21:  Homed  in  Debris  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  57  hours 
prior  to  original  reentry. 
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Figure  4.22:  Summarized  predicted  Debris  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single  A/M 


increase  57  hours  prior  to  original  reentry. 


4-3. 1.2  Second  A/M  Increase.  The  second  point  of  A/M  ratio  increase 
for  the  debris  element  set  began  24  hours  before  the  reference  reentry.  The  purpose 
is  to  be  closer  to  reentry  in  order  to  look  for  a  precision  gain  in  the  time  and  location 
of  the  reentry  prediction. 

Figure  4.23  shows  the  reference  orbit  altitude  of  the  object  as  it  heads  for  reentry 
along  with  the  rest  of  the  predictions.  The  same  behavior  from  the  previous  simulation 
is  observed,  that  is,  all  ten  projections  with  the  A/M  ratio  increase  are  reentered  much 
quicker  and  the  predictions  are  even  closer  together.  The  time  to  reentry  drops  from 
24  hours  to  about  3  hours. 

In  fact,  the  reentry  time  spread  has  dropped  to  45  minutes.  Through  a  closer 
look  at  the  reentry  time  spread,  there  seems  to  be  a  point  that  is  enough  separated 
from  the  group  to  be  considered  an  outlier.  By  removing  this  outlier  the  reentry 
time  spread  lowers  to  16  minutes.  This  can  clearly  be  seen  in  Figure  4.24  where  the 
prediction  using  Atmospheric  Model  Number  1  is  about  half  an  hour  before  the  rest 
of  the  predictions. 

The  drop  in  time  variability  drives  to  a  much  better  reentry  location  spread. 
Figure  4.25  shows  a  plot  of  the  latitude  and  longitude  of  each  reentry  point  and  it  is 
clear  that  the  path  of  reentry  is  much  more  desirable  than  before. 
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Figure  4.23:  Debris  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  24  hours  prior  to  original 
reentry. 
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Figure  4.24:  Homed  in  Debris  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  24  hours 
prior  to  original  reentry. 
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Figure  4.25:  Summarized  predicted  Debris  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single  A/M 


increase  24  hours  prior  to  original  reentry. 


4-3. 1.3  Third  A/M  Increase.  The  last  increase  in  A/M  ratio  was 
chosen  to  be  10  hours  before  reentry.  It  is  expected  to  have  the  best  effect  out  of  all 
previous  points  since  it  is  the  closest  one  to  the  actual  reentry. 

The  same  behavior  from  the  previous  plots  is  observed  in  this  simulation.  The 
simulations  with  the  A/M  increase  dramatically  drop  in  altitude  as  shown  in  Fig¬ 
ure  4.26.  The  time  until  reentry  is  decreased  from  10  hours  to  about  1  hour  and  the 
reentry  time  spread  is  now  down  to  about  20  minutes.  Once  more,  if  we  take  out 
the  first  reentry  prediction  out  of  the  equation,  the  spread  comes  down  to  about  12 
minutes. 

Figure  4.28  shows  a  tightening  of  the  reentry  location  spread  and  the  point  that 
seems  to  be  an  outlier  is  apparent.  Taking  out  this  point,  which  is  the  prediction  by 
Atmospheric  Model  Number  1,  makes  the  reentry  path  cover  36  degrees  in  latitude 
and  just  40  degrees  in  longitude. 

Once  all  3  A/M  increase  simulations  were  performed,  a  plot  displaying  refer¬ 
ence  time  left  to  reentry  as  a  function  of  the  time  prediction  spreads  was  produced. 
Subsequently,  Excel  was  used  to  generate  a  best  fit  trendline  along  with  its  respective 
formula.  The  trendline  equation  has  the  units  of  days  since  dates  were  used  for  all 
the  calculations  in  the  simulations. 

These  first  set  of  simulations  were  propagated  with  the  debris  element  set.  These 
simulations  showed  a  glimpse  of  a  trend.  The  trend  seems  to  be  that  a  quadratic 
relationship  exists  between  the  time  left  to  reentry  and  the  predicted  time  of  reentry 
spread. 


58 


400 


Area/Mass  Ratio  from  0.1  to  1.0 
10  Hrs  Before  Reentry 
vs 

Baseline  Reentry 
(Debris) 


300 


Figure  4.26:  Debris  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  10  hours  prior  to  original 
reentry. 
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Figure  4.27:  Homed  in  Debris  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  10  hours 
prior  to  original  reentry. 
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Figure  4.28:  Summarized  predicted  Debris  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single  A/M 
increase  10  hours  prior  to  original  reentry. 


Reentry  Spread  with  A/M  Increase  from  0.1  to  1.0 
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Figure  4.29:  Summarized  Debris  orbit  orbit  reentry  time  spreads  with  a  single  A/M 
increase  at  57,  24,  and  10  hours  prior  to  original  reentry. 
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4-3.2  ISS  Element  Variability  Change.  The  second  set  of  simulations  for 
this  study  used  the  ISS  orbit  element  set.  The  first  point  of  A/M  ratio  increase  was 
set  7.5  days  prior  to  reference  reentry.  The  change  in  A/M  dropped  the  time  left  to 
reentry  to  less  than  19  hours  but,  had  a  time  reentry  spread  of  over  7  hours.  The 
reentry  location  spread  was  useless  since  the  time  prediction  spread  would  put  it  just 
about  anywhere  in  the  world.  Still,  7  of  the  predictions  were  within  a  trendline. 

The  next  A/M  ratio  increase  was  set  for  24  hours  before  reference  reentry, 
decreasing  the  time  to  reentry  to  less  than  3  hours.  As  in  the  previous  subsection 
this  one  dropped  the  time  reentry  spread  to  43  minutes.  Also,  like  the  previous  set 
of  simulations,  the  spread  seems  to  have  an  outlier  prediction.  Subtracting  it,  gives  a 
reentry  time  spread  of  12  minutes.  Once  more,  the  location  spread  is  very  manageable, 
not  counting  the  outlier  point  which  is  half  an  hour  before  the  rest  of  the  9  points. 

The  last  A/M  ratio  increase  was  set  10  hours  prior  to  the  reference  reentry.  As 
expected,  this  one  gave  the  lowest  reentry  time  spread  at  14  minutes.  When  omitting 
the  prediction  by  Atmospheric  Model  Number  1,  the  time  drops  down  to  7  minutes. 
Also,  the  reentry  location  path  spreads  over  17  degrees  in  latitude  and  28  degrees  in 
longitude.  To  see  further  details  mentioned  above  please  see  Figures  C.l  through  C.9 
in  Appendix  C. 

The  same  quadratic  relationship  that  was  observed  in  the  debris  simulations  on 
Figure  4.29  between  the  time  of  reentry  spread  and  the  reference  time  left  to  reentry 
is  observed  in  the  plot  on  Figure  4.30: 
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Reentry  Spread  with  A/M  Increase  from  0.1  to  1.0 
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Figure  4.30:  Summarized  ISS  orbit  reentry  time  spreads  with  a  single  A/M  increase 
at  180,  24,  and  10  hours  prior  to  original  reentry. 
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4-3.3  Terrasar-X  Element  Variability  Change.  The  last  set  of  simulations 
for  this  section  was  performed  using  an  element  set  from  the  Terrasar-X  satellite. 
This  element  set  not  only  has  a  circular  orbit  but  it  is  also  retrograde.  The  first 
point  of  A/M  ratio  increase  was  set  3  days  before  the  reference  reentry  and  decreased 
the  time  left  to  reentry  to  about  5.5  hours.  It  had  a  reentry  spread  of  over  3  hours. 
This  prediction  spread  was  very  poor  on  the  location  of  reentry  since  only  4  of  the 
predictions  fell  on  a  trendline. 

The  following  point  of  A/M  increase  was  set  24  hours  prior  to  the  reference 
reentry.  The  reentry  time  spread  was  lowered  to  less  than  3  hours.  The  reentry  time 
spread  also  decreased  to  48  minutes.  There  was  only  1  reentry  point  that  was  not 
on  the  trendline.  It  was  half  an  hour  off  the  rest  of  the  predictions.  When  taking 
this  point  off,  the  reentry  time  drops  to  27  minutes  and  the  reentry  locations  have  a 
spread  of  half  a  revolution. 

The  last  point  of  A/M  increase  was  10  hours  prior  to  reentry.  The  time  left  to 
reentry  was  shrunk  down  to  45  minutes.  The  actual  reentry  time  spread  was  of  less 
than  13  minutes.  Finally,  the  reentry  location  path  covered  51.3  degrees  in  latitude 
and  only  10.1  degrees  in  longitude  due  to  its  high  inclination.  One  last  time,  the 
time  reentry  spreads  as  a  function  of  time  left  to  reference  reentry  showed  a  quadratic 
relationship  that  Excel  solved  with  units  of  days.  Further  details  are  in  Figures  C.10 
through  C.18  in  Appendix  C. 

4-3-4  Conclusion  of  Reentry  Variability  Due  to  A/M  Change.  This  study 
established  that  A/M  ratio  manipulation  will  dramatically  change  the  time  left  to 
reentry.  In  addition,  a  great  increase  in  precision  of  the  predicted  reentry  is  obtained 
when  the  manipulation  is  set  closer  to  the  natural  reentry  time.  This  would  be 
done  for  objects  where  predicting  the  place  and  time  of  reentry  are  critical.  Finally, 
the  relationship  between  time  reentry  spread  and  time  left  to  reentry  seemed  to  be 
quadratic.  Upon  a  closer  look  the  relationship  is  actually  much  more  linear  than 
quadratic. 
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Reentry  Spread  with  A/M  Increase  from  0.1  to  1.0 
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Figure  4.31: 


Summarized  Terrasar-X  orbit  reentry  time  spreads  with  a  single  A/M 
increase  at  72,  24,  and  10  hours  prior  to  original  reentry. 
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4-4  Reentry  Prediction  Change  by  A/M  Manipulation  at  a  Single  Point 

The  purpose  of  this  study  is  to  specifically  explore  how  the  reentry  time  predic¬ 
tion  changes  as  we  hold  hold  steady  the  atmospheric  model  and  the  point  where  the 
A/M  manipulation  happens.  This  section  only  uses  the  debris  element  set  to  do  the 
study. 

Figure  4.32  displays  the  reentry  paths  with  A/M  ratios  between  0.01  to  1.0.  It 
can  be  clearly  seen  that  as  you  increase  the  A/M  ratio  the  time  left  on  orbit  decreases. 
For  a  closer  look  at  an  expanded  view  of  the  reentries  with  A/M  ratios  between  0.1 
and  1.0  please  see  Appendix  C,  Figure  C.19.  Also,  Figure  C.20  shows  what  all  reentry 
paths  look  like  compared  with  an  object  that  has  an  A/M  of  0.001.  Excel  is  limited 
to  an  amount  of  points  that  it  is  capable  of  plotting,  therefore,  only  the  last  10  days 
of  data  were  used  in  the  path  of  the  object  with  A/M  of  0.001. 

By  observing  the  plots  in  Figure  4.33,  it  is  clear  there  is  an  exponential  relation¬ 
ship  between  the  A/M  ratio  and  the  time  left  to  reentry.  This  seems  to  stem  from  the 
fact  that  the  atmospheric  density  is  a  factor  in  the  acceleration  equation  3.1.  The  re¬ 
lationship  is  is  better  depicted  by  the  comparison  of  the  time  to  reentry  as  a  function 
of  A/M  plot  against  the  neutral  atmospheric  density  plto  shown  in  Figure  4.34. 

This  is  important,  it  shows  that  for  an  order  of  magnitude  increase  in  A/M  ratio 
an  order  of  magnitude  drop  in  time  left  to  reentry  is  obtained.  Hence,  for  a  Delta  II 
rocket  body  with  a  length  of  6  meters,  a  diameter  of  2.4  meters,  and  a  gross  mass  of 
950  kilograms,  the  maximum  A/M  ratio  possible  would  be  0.015.  Therefore,  if  this 
object  was  in  an  orbit  similar  to  the  one  from  the  USA  193  debris,  instead  of  lasting 
~45  days  in  orbit,  it  would  take  a  balloon  parachute  of  ~  13.5  meters  in  diameter  to 
make  the  rocket  body  reenter  within  5  days. 
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Figure  4.32:  Reentry  predictions  as  A/M  changes  between  0.01  and  1.00  at  a  single  point  in  orbit  path. 
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Figure  4.33:  Summarized  reentry  prediction  times  as  A/M  changes  between  0.001  and 
1.000  at  a  single  point  in  orbit  path. 
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(a)  Time  to  reentry  as  a  function  of  A/M. 


12  14  16  18  20  22  24  26 

Log10  Density  (m  3) 


(b)  Neutral  Atmosphere  (density).  [30] 


Figure  4.34:  Comparison:  Time  to  reentry  as  a  function  of  A/M  vs  Neutral  Density  profile. 


V.  Conclusions  and  Recomendations 


Here  are  the  summarized  conclusions  of  the  studies  as  well  as  proposed  schemes 
for  a  solution.  The  main  proposed  idea  is  to  produce  a  system  that  will  increase 
the  objects  cross-sectional  area.  Consequently,  this  area  increase  would  raise  the 
object’s  drag  and  decrease  its  orbital  lifetime  along  with  gaining  control  over  the 
general  area  where  it  lands. 

5.1  Atmospheric  Unpredictability 

The  first  study  confirmed  how  the  atmospheric  instability  affected  any  predicted 
reentry  that  was  more  than  a  few  days  away  as  it  is  known  in  practice  [38].  It  also 
brought  to  the  front  the  need  to  better  our  atmospheric  modeling  capabilities  when 
it  comes  to  modeling  objects  in  LEO.  These  capabilities  have  increased  in  quality 
with  the  new  NRLMSISE  00,  but  there  is  still  much  room  for  improvement  as  we 
continually  learn  more  about  our  atmosphere’s  changing  nature. 

5.2  Drag  Increase  by  a  Factor 

The  second  study  proved  that  the  increase  in  drag  has  the  effect  of  a  faster 
decay.  The  study  wasn’t  able  to  conclude  whether  there  is  any  difference  in  the  effect 
between  drag  increase  at  apogee  or  perigee.  It  was  observed  that  the  ratio  between 
the  new  time  to  reentry  over  the  old  time  to  reentry  was  inversely  proportional  to  the 
scale  increase  in  drag. 

5.3  Decrease  Variability  in  Reentry  Prediction  and  Manipulation  of 
A/M 

The  third  study  demonstrated  a  quadratic  relationship  between  the  time  left 
to  reentry  and  the  predicted  time  of  reentry  spread,  but  upon  further  scrutiny  the 
behavior  is  basically  linear.  Moreover,  the  third  and  fourth  studies  corroborated  that 
the  increase  in  the  A/M  will  decrease  the  orbital  lifetime  of  the  object.  This  will 
make  it  possible  for  objects  to  reenter  sooner.  It  is  important  to  remember  all  studies 
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were  done  using  orbital  element  sets  that  have  altitudes  not  much  greater  than  500 
kilometers  (the  lower  end  of  LEO  orbits).  This  was  set  with  the  purpose  that  any 
object  that  could  use  the  system  being  proposed  would  have  a  periapsees  low  enough 
to  the  earth  where  it  can  make  use  of  the  “thicker”  atmospheric  density. 

Not  only  can  a  system  that  increases  the  object’s  A/M  be  brought  down  sooner, 
it  will  consequently  gain  precision  in  the  predicted  reentry  time.  It  was  confirmed  that 
the  reentry  prediction  decreases  greatly  in  variability  as  it  comes  closer  to  the  actual 
reentry.  Therefore,  when  forcing  an  object  to  reenter,  you  automatically  gain  precision 
in  the  predicted  reentry  time.  Along  with  variability  reduction  in  time  of  reentry,  the 
next  benefit  is  an  actual  decrease  in  variability  of  the  predicted  reentry  location. 

All  simulations  decreased  their  variability  in  location  by  setting  their  A/M  in¬ 
crease  closer  to  reference  reentry.  Even  the  prediction  with  the  A/M  increase  one  day 
out  gave  about  a  one  revolution  variability,  or  what  it  amounts  to  a  “ring  possibility” 
of  reentering  around  the  earth.  This  ring  maybe  enough  of  a  prediction  for  the  reen¬ 
tering  object  if  it  is  known  to  have  minimal  chances  of  survival.  It  was  also  shown 
that  even  more  precision  may  be  gained  if  the  A/M  increase  is  postponed  to  be  closer 
to  actual  reentry. 

For  this  reason,  it  is  also  suggested  the  development  of  software  that  can  solve 
for  a  needed  cross-sectional  area  (or  balloon  size)  given  a  desired  window  of  reentry. 
The  software  would  need  inputs  such  as  the  object’s  physical  properties  and  its  orbital 
parameters.  It  would  require  a  robust  dynamic  atmospheric  model.  The  inspiration 
for  this  idea  came  from  the  fourth  study,  where  it  was  concluded  that  the  time  to 
reentry  has  a  relationship  that  is  exponential  to  changes  in  A/M.  It  was  determined 
that  this  behavior  was  due  to  the  density  of  the  Neutral  Atmosphere.  This  is  because 
atmospheric  density  is  a  factor  in  the  object’s  acceleration  due  to  drag. 


72 


5.4  Parachute  Proposal 

The  proposed  approach  to  increasing  the  A/M  is  through  the  on-orbit  inflation 
of  a  balloon.  The  orbiting  object  will  carry  an  attached  package  that,  when  signaled, 
will  activate  and  simply  open  a  compartment  where  a  bag  will  be  exposed.  This  bag 
will  be  deployed  as  it  is  inflated  and  will  start  acting  as  a  parachute.  Another  idea 
for  increasing  the  A/M  would  be  through  the  deployment  of  a  sail.  The  problem  with 
this  idea  is  that  it  would  probably  need  a  more  sophisticated  system  for  deployment 
and  require  an  attitude  control  system. 

5-4-1  Design  Parameters.  It  is  desirable  to  make  the  apparatus  as  cheap 
and  simple  as  possible.  It  is  well  known  that  the  more  complicated  an  object  becomes, 
the  more  things  that  can  go  wrong  with  it.  Therefore,  from  the  inception  of  the  idea, 
simplicity  in  design  must  be  aggressively  followed. 

A  minimum  amount  of  functions  in  the  apparatus  is  a  must  while  minimizing 
moving  parts.  The  single  most  complicated  idea  at  this  point  would  be  the  desire  to 
have  its  own  power  supply.  The  goal  is  for  the  system  to  have  autonomy  and  be  able 
to  come  into  action  when  needed.  If  it  has  its  own  power  source,  it  may  be  able  to 
force  the  reentry  of  a  malfunctioned  satellite. 

Keeping  the  design  simple  will  add  robustness  and  drive  to  a  low-mass  end- 
product.  This  will  keep  the  cost  down  for  launch,  while  a  separate  power  source  will 
provide  the  system  with  versatility  to  go  on  different  orbiting  objects.  The  list  would 
include  satellites,  rocket-bodies,  platforms,  and  any  big  enough  pieces  of  debris  where 
it  may  be  used. 

5-4-2  Prototype  Use.  A  satellite  parachute  prototype  would  probably  be 
better  used  on  a  fuel  tank  which  will  already  be  left  in  orbit  such  as  the  one  from  a 
Delta  II  rocket  body  second  stage.  It  may  also  be  used  on  experimental  satellites. 

Then  again,  the  use  of  the  prototype  on  a  more  simple  object  such  as  a  rocket 
body  could  constitute  faster  deployment.  The  reason  for  this  is  that  its  power  source 


73 


may  not  be  as  complicated  as  one  that  has  to  make  sure  the  prototype  survives  until 
the  end  of  a  satellite’s  mission.  It  does  mean  that  some  modifications  must  be  made 
to  the  launch  system  (rocket  body).  It  must  be  kept  in  mind  that  this  is  to  further 
the  cause  of  increasing  our  knowledge,  experience,  and  practical  expertise  of  a  new 
system  that  will  help  alleviate  the  debris  problem  in  space  and  further  protect  the 
precious  resource  that  is  our  limited  orbital  space. 

5-4-3  Further  Research  Needed.  Pageos  I,  Echo  I,  and  Echo  II  satellites 
proved  that  balloons  are  capable  of  surviving  for  long  periods  of  time  up  in  space 
even  when  they  are  intercepted  by  micrometeors.  Still,  a  small  collision  probability 
assessment  must  be  performed  since  there  is  a  considerable  larger  amount  of  manmade 
debris  in  space. 

Many  of  the  reports  from  the  missions  of  the  balloon  satellites  may  be  exploited 
to  benefit  this  project.  It  is  important  to  remember  that  these  systems  are  expected 
to  go  through  a  higher  level  of  stress  than  previous  systems  did.  The  reason  for  this 
is  because  their  capability  to  decrease  the  orbital  lifetime  of  objects  will  be  tested. 

As  a  final  note,  it  is  important  to  realize  that  the  simulations  only  varied  the 
A/M  and  kept  the  coefficient  of  drag  constant.  This  is  not  entirely  accurate.  When 
increasing  the  geometry  of  the  object,  the  coefficient  of  drag  is  also  affected.  Therefore, 
the  effect  of  the  increase  in  surface  area  should  actually  be  greater. 
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Appendix  A.  Orbital  Element  Sets 


A.l  Orbital  Element  Sets  Used  For  Simulations 

The  following  figures  provide  the  orbital  element  sets  that  were  used  as  the 
beginning  points  for  the  simulations. 

They  give  a  quick  description  and  international  designator  of  the  orbiting  object 
from  which  they  were  acquired  along  with  the  web  address  where  the  elements  were 
obtained. 

The  last  figure  in  this  Appendix  is  a  set  of  keys  to  decipher  element  sets.  The 
bottom  key  is  the  one  to  be  used  for  the  element  sets  presented. 
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International  Designator  2006-057 


2006-057A 

USA  193  is  an  American  military  satellite  that  was  launched  from  Vandenberg  AFB  at  21:00  UT  on 
2006  December  14.  It  is  a  highly  classified  spacecraft,  owned  and  operated  by  the  National 
Reconnaissance  Office  (NRO).  No  further  details  are  available. 

Debris  SCC  number  29651 

The  satellite  was  destroyed  on  Feb  21, 2008. 

Accessed  on  23  July  2008  at  1647  ET 
5?  From 

http://celestrak.com/NORAD/elements/usa-193-debns.txt 

USA  193  DEB 

1  32601D  06057BY  08204 . 81D39S 92  .00270106  32405-4  65292-3  0  2673 

2  32601  56.5124  127.6791  C163670  250.1491  108.2035  15.75949212  22370 


Figure  A.l:  Orbital  element  set  for  Debris  orbit. 


International  Designator  1998-067 

1998-067A 

ISS-ZARYA  is  a  Russian-built  first  module  of  the  International  Space  Station  (ISS)  that  was 
launched  by  a  Proton-K  rocket  from  Baikonur  at  06:40  UT.  The  ISS  is  also  known  by  the  Russian 
acronym,  MKS  standing  for  the  transliterated  Mezhdunarodnii  Kosmicheskii  Stantsii.  It  will  require 
43  rocket  launches,  including  about  35  Shuttles,  over  at  least  a  five  year  period  to  complete  the  16- 
country,  60  m  x  24  m  x  21 .5  m,  454-ton,  1 1 0-kW,  and  >$50  billion  ISS  (or  whatever  its  eventual 
name/acronym  may  turn  out  to  be).  ZARYA  ("Dawn")  module  has  a  mass  of  27  tons  and  will 
contribute  power,  attitude  control,  fuel,  and  command/control  coordination  to  the  other  33  (USA  18, 
Russia  9,  Japan  3,  ESA  2,  and  Canada  1)  modules.  No  definitive  information  is  currently  available 
regarding  ISS's  scientific,  commercial,  or  military  vision.  It  is  likely  that  some  or  all  of  the  science 
payload  now  housed  in  the  140-ton  Mir  station  may  be  transferred  to  the  ISS  in  the  event  of  having 
to  dissemble  and  deorbit  Mir;  the  orbital  parameters  of  the  two  are  similar.  Initial  orbital  parameters 
of  ISS-ZARYA  were  period  92  min,  apogee  396  km,  perigee  384  km,  and  inclination  51.6  deg. 

SCC  number  25544 

Accesed  on  30  Jul  08  at  1234  ET 

From  http://celestrak.com/NORAD/elements/stations.txt 

ZARYA 

1  25544U  98067A  08211.87215691  .00006360  00000-0  50579-4  0  9514 

2  25544  51.6417  156.9681  0010004  11.4376  39.1422  15.72722297  555256 


Figure  A. 2:  Orbital  element  set  for  ISS  orbit. 


2007-026A 


International  Designator  2007-026 


TERRASAR-X  is  a  German  (DLR)  Synthetic  Aperture  Radar  that  was  launched  by  a  Dnepr  rocket 
from  Baikonur  at  02: 1 4  UT  on  2007  June  1 5.  It  will  map  the  Earth  and  oceans  at  a  rate  of  one 
million  square-km/day,  provide  scientific  data  such  as  sea  ice  coverage,  vegetation/crop  estimates, 
and  military  reconnaissance  information,  all  at  1.0-m  resolution. 

The  initial  orbital  parameters  were  period  95  min,  apogee  510  km,  perigee  507  km,  and  inclination 
97.5°. 


Accessed  on  2  December  2008  at  13:36  ET 
From 

http://celestrak.eom/satcat/2007/2007-026.asp#A 

http://celestrak.com/NORAD/elements/resource.txt 

TERRASAR-X 


1  31698U 


2  31698 


07026A 


97 . 4442 


08337.19636624 


341 . 4278 


.00000629 


0001424 


00000-0 


77.7037 


33137-4  0  6377 


52.2730  15.19157484  8136 


Figure  A. 3:  Orbital  element  set  Terrasar-X  orbit. 
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Figure  A. 4:  Orbital  element  set  keys. 


Appendix  B.  STK  Propagator  Inputs 

B.l  Set  Up  Inputs  for  Simulations 

The  following  figures  are  the  propagator  inputs  used  for  the  initial  runs  of  the 
simulations  of  all  three  orbits. 

It  is  important  to  remember  that  the  inputs  changed  as  different  simulations 
were  required  by  the  different  studies  or  to  have  the  propagator  produce  better  results. 

Such  changes  included  A/M  ratio  manipulations,  integrator  time  step  manipu¬ 
lation  and  atmospheric  model  type  changes. 
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Educational  Use  Only 
Satellite-MSISE  1990 

Propagator  Initial  Conditions 


Propagator  Name  =  HPOP 

Orbit  Epoch  =  22  Jul  2008  19:26:00.000000000 

Position  -  J2000  -  x  =  -4138831.51430963 
Position  -  J2000  -  y  =  5341288.19322933 
Position  -  J2000  -  z  =  17736.34682893 
Velocity  -  J2000  -  x  =  -3226.12586528 
Velocity  -  J2000  -  y  =  -2367.10813792 
Velocity  -  J2000  -  z  =  6530.86213070 

Mass  =  1000.000 

Gravity  File  Name  =  C:\Program  Files\AGI\STK  8\STKData\CentralBodies\Earth\WGS84_EGM96.grv 
MaxDegree  =  50  MaxOrder  =  50 
Include  Solid  Tides?  No 
Include  Ocean  Tides?  No 

Include  Drag?  Yes 
Drag  Coefficient  =  2.0000 
Area/Mass  Ratio  =  0.1000 
Atmos  density  model  =  MSISE  1990 
Use  Flux  File?  Yes 

Flux  Filename  =  C:\Program  Files\AGI\STK  8\DynamicEarthData\SpaceWeather-vL2.txt 
Use  Approximate  Altitude  for  Drag?  Yes 
Use  Apparent  Sun  for  Drag?  Yes 

Include  Solar  Radiation  Pressure?  No 
Include  general  relativity  correction?  No 
Include  albedo?  No 
Include  thermal  radiation  pressure?  No 

Integrator  =  Runge-Kutta-Fehlberg  7-8 

Relative  Error  Tolerance  =  1.000000e-013 

Integrator  Step  Size  =  Variable 

Report  On  Fixed  Step?  Yes 

Equations  of  Motion  =  VOP  Universal  Variables 

Independent  Variable  =  Time 

Interpolation  Method  =  Lagrangian  VOP  (mu  =  3.98600441 80000000e+0 14) 

Interpolation  Order  =  8 

Integrator  stops  when  altitude  is  below  0.0000  m 
Computing  partials?  No 


Figure  B.l:  STK  propagator  inputs  for  Debris  orbit  simulation. 
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Educational  Use  Only 
Satellite -MSIS  1990 

Propagator  Initial  Conditions 


Propagator  Name  =  HPOP 

Orbit  Epoch  =  1  Aug  2008  15:21:00.000000000 

Position  -  J2000  -  x  =  -5193760.32542793 
Position  -  J2000  -  y  =  -1297506.48023412 
Position  -  J2000  -  z  =  4076421.27492316 
Velocity  -  J2000  -  x  =  4285.68711 167 
Velocity  -  J2000  -  y  =  -51 17.57790620 
Velocity  -  J2000  -  z  =  3832.28596309 

Mass  =  1000.000 

Gravity  File  Name  =  C:\Program  Files\AGI\STK  8\STKData\CentralBodies\Earth\WGS84_EGM96.grv 
MaxDegree  =  50  MaxOrder  =  50 
Include  Solid  Tides?  No 
Include  Ocean  Tides?  No 

Include  Drag?  Yes 
Drag  Coefficient  =  2.2000 
Area/Mass  Ratio  =  0.1000 
Atmos  density  model  =  MSISE  1990 
Use  Flux  File?  Yes 

Flux  Filename  =  C:\Program  Files\AGI\STK  8\DynamicEarthData\SpaceWeather-vL2.txt 
Use  Approximate  Altitude  for  Drag?  Yes 
Use  Apparent  Sun  for  Drag?  Yes 

Include  Solar  Radiation  Pressure?  No 
Include  general  relativity  correction?  No 
Include  albedo?  No 
Include  thermal  radiation  pressure?  No 

Integrator  =  Runge-Kutta-Fehlberg  7-8 
Relative  Error  Tolerance  =  1.000000e-013 
Integrator  Step  Size  =  Variable 
Report  On  Fixed  Step?  Yes 
Equations  of  Motion  =  Cowell 
Independent  Variable  =  Time 
Interpolation  Method  =  Lagrangian 
Interpolation  Order  =  7 

Integrator  stops  when  altitude  is  below  0.000  m 
Computing  partials?  No 


Figure  B.2:  STK  propagator  inputs  for  ISS  orbit  simulation. 
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Educational  Use  Only 
Satellite-Satellite  1 

Propagator  Initial  Conditions 


Propagator  Name  =  HPOP 

Orbit  Epoch  =  2  Dec  2008  04:42:46.043136000 

Position  -  J2000  -  x  =  -4412479.14888076 
Position  -  J2000  -  y  =  761502.76035328 
Position  -  J2000  -  z  =  523 1259.08903455 
Velocity  -  J2000  -  x  =  -5324.63669209 
Velocity  -  J2000  -  y  =  2457.28610748 
Velocity  -  J2000  -  z  =  -4847.81515618 

Mass  =  1000.000 

Gravity  File  Name  =  C:\Program  Files\AGI\STK  8\STKData\CentralBodies\Earth\WGS84_EGM96.grv 
MaxDegree  =  50  MaxOrder  =  50 
Include  Solid  Tides?  No 
Include  Ocean  Tides?  No 

Include  Drag?  Yes 
Drag  Coefficient  =  2.0000 
Area/Mass  Ratio  =  0.1000 
Atmos  density  model  =  MSISE  1990 
Use  Flux  File?  Yes 

Flux  Filename  =  C:\Program  Files\AGI\STK  8\DynamicEarthData\SpaceWeather-vL2.txt 
Use  Approximate  Altitude  for  Drag?  Yes 
Use  Apparent  Sun  for  Drag?  Yes 

Include  Solar  Radiation  Pressure?  No 
Include  general  relativity  correction?  No 
Include  albedo?  No 
Include  thermal  radiation  pressure?  No 

Integrator  =  Runge-Kutta-Fehlberg  7-8 

Relative  Error  Tolerance  =  1.000000e-013 

Integrator  Step  Size  =  Variable 

Report  On  Fixed  Step?  No 

Equations  of  Motion  =  VOP  Universal  Variables 

Independent  Variable  =  Time 

Interpolation  Method  =  Lagrangian  VOP  (mu  =  3.9860044180000000e+014) 

Interpolation  Order  =  8 

Integrator  stops  when  altitude  is  below  0.0000  m 
Computing  partials?  No 


Figure  B.3:  STK  propagator  inputs  for  Terrasar-X  orbit  simulation. 
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Appendix  C.  Graphs  From  Chapter  Four 


C.l  Graphs  From  Section  Three  and  Four  of  Chapter  Four 

Since  the  amount  of  figures  in  sections  3  and  4  of  Chapter  4  were  impeding  a 
better  reading  flow  for  the  thesis,  they  were  moved  to  this  appendix. 

The  graphs,  along  with  their  summaries,  depict  the  predicted  reentry  paths  of 
the  different  orbits. 
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Figure  C.l:  ISS  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  180  hours  prior  to  original 
reentry. 
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Figure  C.2:  Homed  in  ISS  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  180  hours  prior 
to  original  reentry. 


A/M  Scale  Increase  Point 
_ 180  hrs _ 

Reentry  Time  Spread 
(hrs:mm:ss) 
7:11:00 


Sorted  by  Date 


Atmospheric 

Model 

Reentry  Time 
Time 

Latitude 

(deg) 

Longitude 

(deg) 

Altitude 

(km) 

1 

8/10/08  0:15 

37.31 

294.30 

106.42 

4 

8/10/08  2:22 

-51.07 

41.53 

104.65 

3 

8/10/08  6:28 

2.36 

234.51 

97.07 

6 

8/10/08  6:52 

-51.32 

-30.05 

99.95 

9 

8/10/08  6:54 

-50.52 

-23.44 

98.32 

10 

8/10/08  6:58 

-44.47 

-1.99 

100.41 

8 

8/10/08  6:59 

-42.58 

2.05 

98.13 

7 

8/10/08  7:06 

-21.82 

28.12 

99.16 

5 

8/10/08  7:20 

18.28 

58.24 

99.71 

2 

8/10/08  7:26 

36.18 

76.83 

101.45 

oo 

-n 


Reentry  Locations 

A  ▲ 

Latitude 

(deg)  Q 

/  ; ;  A 

-35 

15  /  65  115  165  215  265  315 

/  Longitude  (deg) 

A 

- 89— 

Figure  C.3:  Summarized  predicted  ISS  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single  A/M 


increase  180  hours  prior  to  original  reentry. 
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Figure  C.4:  ISS  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  24  hours  prior  to  original 
reentry. 
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Figure  C.5:  Homed  in  ISS  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  24  hours  prior 
to  original  reentry. 
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Figure  C.6:  Summarized  predicted  ISS  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single  A/M 


increase  24  hours  prior  to  original  reentry. 
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Figure  C.7:  ISS  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  10  hours  prior  to  original 
reentry. 
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Figure  C.8:  Homed  in  ISS  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  10  hours  prior 
to  original  reentry. 
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Figure  C.9:  Summarized  predicted  ISS  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single  A/M 


increase  10  hours  prior  to  original  reentry. 
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Figure  C.10:  Terrasar-X  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  72  hours  prior  to 
original  reentry. 
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Figure  C.ll:  Homed  in  Terrasar-X  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  72  hours 
prior  to  original  reentry. 
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Figure  C.12:  Summarized  predicted  Terrasar-X  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single 


A/M  increase  72  hours  prior  to  original  reentry. 
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Figure  C.13:  Terrasar-X  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  24  hours  prior  to 
original  reentry. 
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Figure  C.14:  Homed  in  Terrasar-X  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  24  hours 
prior  to  original  reentry. 
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Figure  C.15:  Summarized  predicted  Terrasar-X  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single 


A/M  increase  24  hours  prior  to  original  reentry. 
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Figure  C.16:  Terrasar-X  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  10  hours  prior  to 
original  reentry. 
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Figure  C.17:  Homed  in  Terrasar-X  orbit  reentry  predictions  by  ten  atmospheric  models  with  a  single  A/M  increase  10  hours 
prior  to  original  reentry. 
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Figure  C.18:  Summarized  predicted  Terrasar-X  orbit  reentry  times  and  locations  by  ten  atmospheric  models  with  a  single 


A/M  increase  10  hours  prior  to  original  reentry. 
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Figure  C.19:  Homed  in  Reentry  predictions  as  A/M  changes  between  0.10  and  1.00  at  a  single  point  in  orbit  path. 
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Figure  C.20:  All  reentry  predictions  as  A/M  changes  between  0.001  and  1.000  at  a  single  point  in  orbit  path. 
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